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ABSTRACT 
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and analysis. In addition, the internal balance that vas 
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used and programs developed for data acquisition and 
analysis were adequate. However, balance calibration 
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the roll component strain gage. Both of these problems vill 
have to be corrected before accurate readings can be 
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LI.  HTRODUCETON 


A. BACKGROUND 

This project used the 3.5" x 5.0" lov speed vind tunnel 
at the Naval Postgraduate School to continue the 
aerodynamic study of the effects of helicopter noses and 
tails on drag conducted by Major Mair [Ref. 1] and Major 
Sargent (Ref. 21. Their work included the construction of 
nine various helicopter configurations, a sting mounted 
support system and a four-degree-of-freedom balance that 
was internal to the model. They also studied the airflow 
around the various configurations with cotton tufting. 

To reduce and analyze the data, Majors Mair and Sargent 
developed several independent computer programs. These 
programs calibrated the balance, corrected the raw data 
from the wind tunnel runs, computed the equivalent flat 
plate area and produced plots of various parameters for 
comparison. However, due to problems encountered with the 
drag component of the balance, their results were 


inconclusive. 


B. GOALS 

The primary goal of this project was to provide the 
students enrolled in the helicopter design class at the 
Naval Postgraduate School with a laboratory type set-up 


that would develop realistic Equivalent Flat Plate Area 


information for various helicopter configurations. This 
included creating interactive computer programs that the 
students could use with the personal computer located at 
the wind tunnel. 

The secondary goals of this project were to upgrade the 
two-axis internal wind tunnel balance to a six-axis 
internal balance and to develop a calibration rig to 
evaluate the interactions of the balance components. 

Three landing gear configurations were also designed 
and constructed for future analysis of the drag that they 


add to the helicopter. 
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no APEROACH TONTHE. PROBLEM 


A. LANDING GEAR DESIGN 

To provide a realistic representation of landing gear 
used with modern-day helicopters, one type of landing gear 
was selected for each nose Shape (Figures 2.1-2.3). 
Planviews for the three types of landing gear were prepared 
and are included in Appendix A. 

For the attack nose, a skid type of landing gear was 
constructed of aluminum tubing secured to an aluminum 
plate. This type of landing gear is considered a fixed gear 
but was chosen because of its vide use for numerous 
helicopters. For both the smooth nose and blunt nose a 
simulated retractable gear was constructed of wood 
stubvings and model airplane tires. Threaded inserts were 
mounted in the noses and stubwings to allow easy removal of 
the wheel and strut assemblies. This allowed the models to 


be tested in both the clean and dirty configuration. 


B. INTERNAL BALANCE MODIFICATION 

A modified Mair/Sargent balance, Figures 2.4-2.7, was 
used for this project: 

TOW-improveq tlie 7 outbut recorded GEOM the axial 
component, the cuts (see Figure 2.4) alongside that cavity 
were increased by 1/16 of an inch. In addition, the cavity 


itself was squared off thus reducing the curvature of the 


jer 





Figure 2.1 Attack Nose vith gear 





Figure 2.2 Smooth Nose vith gear 





Rigure 2.3 Blunt Nose vith gear 
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Elmin Re-Wired Internal Balance 
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surfaces in the cavity. These modifications increased the 
flat area upon which the strain gages could be placed. 

An area (see Figure 2.6), similar to that for the 
pitching moment component, was cut to allow recording of 
the yawing moment component. 

For compatability and increased accuracy, the aluminum 
strain gages were removed and replaced by EA-09-062AQ-350 
stainless steel ones.- These gages were smaller and thus 
allowed for better placement within the cavities. To 
provide a longer life, the gages were bonded = the balance 
with an M-Bond AE-15 adhesive system. They were cured at a 
temperature of 150 degrees Fahrenheit for two hours. 

It was desired to record the six component forces on 
the helicopter; lift, drag, yaw, pitching moment, yawing 
moment and rolling moment. However, since the internal 
balance was designed to record only four components, the 
sting support was instrumented to record the yaw force and 
rolling moment (Figure 2.8). 

The gages for the yaw component were placed on the side 
of the sting support to undergo tension and compression 
when subjected to a yawing force. The gages for the rolling 
moment component were placed on top and bottom of the sting 
Support at a 45 degree angle to the sting axis. Thus, they 
experienced tension and compression when the model vas 


Subjected to a rolling moment. 
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Figure 2.6 Gages 





Figure 287 Drag Component Gage 





F ici rew2.8 Yaw and Roll Moment Gages 
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Ce BALANCE CALIBRATION 

In order to determine any interaction between the six 
recording components, a calibration rig was designed and 
mounted to the main fuselage (Figure 2.9). The pans and 
pullies were arranged to simulate the twelve forces and 
moments that the helicopter models would experience in the 
wind tunnel. Each pan, or component was loaded from zero to 
twenty pounds in one pound increments. With each loading, 
readings from the other five components were recorded. This 
method was conducted until all six components were loaded 
in both the positive and negative direction and produced 
thirty interaction matrices. 

A balance calibration progran, Figure A.16, was written 
to determine the relations required to convert the raw data 
counts to actual forces and moments. For each component 
loading, the prime gage constants were determined using the 


following least squares curve fit [Ref. 3]: 


[Sum(Xi)*2]*Kl + [Sum(Xi*3) ]*K2 = Sum(Xix*Yi) 
(Sum(Xi)731“K1 + [Sum(Xi*4)]*K2 = Sum((Xi*2)*Yi) . 
where Xi is the raw data count 
and Yi is the applied load 
Once these constants were computed for all twelve loadings, 
the raw data counts were converted to forces and moments. 
Then for each interaction matrix, the same least squares 
equations were used to determine the interaction 


coefficients. These prime gage constants and interaction 
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coefficients vere stored for use in the reduction of the 


raw data recorded during the wind tunnel tests. 


D. DATA ACQUISITION 

To improve data reduction, it was desired to have a 
data acquisition system built into the wind tunnel 
instrumentation. The strain gages on the internal balance 
and support sting were connected to a Pacific Instruments 
model 8255 transducer amplifier (Figure 2.10). The outputs 
from this unit were then routed to data acquisition cards 
mounted in an IBM PC AT. A data acquisition program, Figure 
A.10, was incorporated into the interactive program set 
that allowed for quick and easy recording of test data. 

Having the transducer amplifier connected to the strain 
gages allowed for easy zeroing and balancing of the gages. 
The data acquisition set-up provided real time readout of 
the forces that the helicopter was experiencing in the wind 


tunnel. 
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Data Acquisition Equipment 
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III. SOLUTION 1065155. BEIN 


A. DATA COLLECTION AND REDUCTION 

The data collection program, Figure A.10, vas designed 
to record the data in a standardized collection method. 
Interactive steps insiet the user to follow the 
procedures listed below. 

l). Zero all amplifiers without the model in place. 


2). Mount the model to the sting support and record 
a zero force reading. 


3). Place the calibration switches to the 4 position 
and record a calibration reading. f 

4). Replace the calibration switchres to the center 
position and start the wind tunnel. 

5). Record first and last data point at zero angle of 


attack with all other readings between plus ten 
degrees and minus eight degrees angle of attack. 

By using a delay loop, changes could be made to the 
angle of attack of the model before the next set of data 
points were recorded. For each tunnel speed, the data 
recorded for the various angles of attack vere stored on 
separate files. 

Once the tunnel runs were completed and the data was 
recorded on file, another program, Figure A.11, was 
designed to convert the raw readings to readings of forces 
and monemts. The sixty interaction coefficients developed 
from the balance calibration were incorporated into this 


data reduction. The following equations, obtained from Mr. 
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David Backs at the NASA AMES Reseasrch Center, vere used to 
correct for the interaction of the balance components. 
Fl' = Fl - (delFl/delF2)*F2 - (delFl/delF3)*F3 - .......... 
ー (de1F1/de1F2 へ ^2) *F2^2 = ....... -(delFl/delF6*2) *F6^2 
Pu rı = (ии /ДеЛЕ2)*Е28=- 1... әәә ә ии 
- (delFl/delF202)*F2'42 ....... -(de1F1/de1F6^2 ) *F6! へ 2 
These equation were written for each component and 
placed in an iterative loop that checked the difference 
between the two prime values. After the forces were 
corrected for component interaction, a weight tare equation 


was used to correct for the weight of the model 


B. ANALYSIS OF DATA 

The main file of the interactive programs, Figure A.9, 
controlled the data recording, data reduction and data 
analysis with a menu format. For data analysis, the stored 
converted data readings were used to calculate the lift and 
drag coefficients and equivalent flat plate area for each 
angle of attack. Files were created for coefficient of lift 
versus angle of attack, coefficient of drag versus 
coefficient of lift, coefficient of drag versus coefficient 
of lift squared and equivalent flat plate area versus angle 
of attack. There was one file of each created for each 
tunnel speed. 

A plotting routine, Figure A.15, was included in the 
interactive programs to allow quick analysis of the 


recorded data. To aid in the analysis, up to three plots 


21 


could be shovn on one graph. A delay loop vas incorporated 
into the plotting “routine to ALLO VM SE MM TS OPE FO NO l 


obtaining a hard copy of the graph by using the print 


screen command. 


22 


IV. RESULTS 


A. BALANCE CALIBRATION 

The balance vas loaded in the tvelve component 
directions and the balance calibration program vas used to 
produce the tables B.2 through B.13. During the balance 
calibration, a drift in each of the components vas noted 
under steady state conditions. To correct for this problem 
it was assumed that each 27:50: had a constant drift 
rate. The difference between the first and last zero 
reading was divided by the number of data points taken. 
This correction factor was then applied to each data point. 
This method produced good correlation between different 


(data points taken for the same load. 


B. DATA COLLECTION 

Test runs were conducted te evaluate the data 
acquisition program. The vibrations of the model in the 
vind tunnel resulted in erratic fluctuations in the sense 
Hüidğicators ef the amplifiers. To correct for this problem, 
the acquisition program was modified to collect one hundred 
samples at a rate of five hundred samples per second. The 
average of these one hundred samples was taken as one data 
point. This method produced constant readings for data 


points of similar conditions. 


ә 


C. DATA REDUCTION AND ANALSYIS 

The data reduction program was used to convert the raw 
test data. The interaction equations diverged instead of 
converging to a single value. This pointed out a problem 
vith the interaction coefficients. Examination of the 
calibration tables revealed extremely large interactions 
betveen the loading of the lift component and reaction in 
- the roll moment component. There was also noted a large 
interaction between the loading of the pitch moment 
component and reaction in the lift component. 

The reduction program was modified to correct for the 
interactions between the lift and drag components only. The 
reduction program was again executed using the recorded 
test data. This time the interactions converged. This 
proved the validity of the interaction equations and 
confirmed the problems with the balance calibration. 

A test file was created to display the plotting 
capabilities of the interactive programs. Figure 4.1 shows 


the results. 
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Figure 4. 


V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The data acquisition system and programs developed were 
adequate for recording the forces and moments experienced 
in the wind tunnel. The extremely large interactions 
between the lift and roll moment components and between the 
pitch moment and lift components (see Tables B.2 - В.13) 
precluded any analysis of wind tunnel data. However, the 
data acquisition system and balance calibration program 
warrant further research and development of the balance 
system. 

The interactive programs greatly reduced the workload 
required in the data acquisition and analysis phases. They 
provided a quick and easy means for the analysis of the 


recorded data. 


B. RECOMMENDATIONS 
The following are given as recommendations to improve 


the calibration of the internal balance system. 


1. Balance Modification 
The location of the roll moment strain gage was 
decided to be the cause of the large interaction between 
the lift and roll moment component. This gage was located 


on the back end of the sting support system. The roll 


26 


moment component vas calibrated by applying a torque to the 
center section of the model that vas mounted to the 
internal balance. Vhen the lift component vas loaded, a 
large bending moment vas felt by the roll moment gage and 
thus producing the large interactions. 

By placing the roll moment strain gage at a 
forty-five degree IO the same cut-out section as the 
pitch moment gage, the torque applied to the center section 
can be used for calibration. Also, the loading of the lift 
component will not greatly affect the reaction of the roll 


moment component. 


2. Calibration Rig Modification 
By using the center section of the model to mount 
the calibration pans for the pitch moment component, the 
strains were incorrectly transmitted to the lift component. 
A separate calibration set-up will have to be designed that 
isolates the force appplied to the pitch moment component 


from the rest of the balance. 
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APPENDIX A 


SKETCHES AND” PROGRAMS 





Figure A.1 Nose Gear for Smooth Nose 
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Figure A.2 Main Gear for Smooth Nose 





Figurer 


Stubwing for Smooth Nose (Top View) 
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Skid Gear and Wing for Attack Nose 


Figure A.4 
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Flgure A.5 Skid Gear for Attack Nose (Top Viev) 





Figure A.6 Nose Gear for Blunt Nose 
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Main Gear for Blunt Nose 


Figure A.7 
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Figure A.8 Stubving for Blunt Nose (Top Viev) 
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KEM=== mz MAIN.BAS (MAIN CONTROL PROGRAM)----- 
: PATRICK A. VITT 
I 20 JULY r225 


FOR I = 1I TO 2 
= 1 THEN FILL 31 ELSE "FZZ 
IF I = 1 THEN EDGE 9 ELSE EDGE = 10 
= 1 THEN X1=106: X2=0: X3=56 ELSE X1=108: X2=2: 


IF I = 1 THEN Y1=48: Y2=94: Y3=140 ELSE Y1=46: Y2=92: 
Y3=138 

' ир" 

DRAW "C=EDGE; BM=X1; ,=Y1;U30R4D12R12U12R4D30L4U12" 
DRAW "L12D12L4BE1P=FILL; ,=EDGE;'" 

' "Z" 

X1= X1 + 26 

DRAW ""C=EDGE;BM=X1; ,=Y1;U30R20D4L16D9R8D4L8D9R16D4L2 
OBE1P=FILL; ,=EDGE;" 

' "TW 

X1= X1 + 26 

DRAW "C=EDGE; BM=X1;,=Y1;U30R4D26R16D4L20BE1P=FILL;, 
=EDGE;" 

' "o" 

Xl- XI + 26 

DRAW "C=EDGE; BM=X1; ,=Y1;U30R20D30L20BR4BU4U22R12D22L 
12BG1P=FILL; ,=EDGE;" 

'NEXT LINE 

' m. 

DRAW "C=EDGE; BM=X2; ,=Y2;U30R20D4L16D22R16D4L20BE1P= 
FILL; ,=EDGE;" 

' "o" 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2; ,=Y2;U3O0R20D30L20BR4BU4U22R12D22L 
12BG1P=FILL; ,=EDGE;" 
, UN! 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2;,=Y2;U30R4F12U12R4D30L4H12D12L4BE 
1P=FILL;,=EDGE;" 

' "p" 

X2= X2 + 26 

DRAW "C=EDGE ;BM=X2; ,=Y2;U30R20D4L16D6R8D4L8D16L4BE1P= 
FİLE? SEDGE ? " 

' "(T 

X2= X2 + 26 

DRAW "C=EDGE; BM=X2; ,=Y2;U30R4D30L4BE1P=FILL; ,=EDGE;" 
1 Med 


Figure A. 9 MAIN.BAS - Main Controlling Program 
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770 
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840 
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870 
880 
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900 
210 
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X2= X2 + 14 

DRAW "C=EDGE; BM=X2; ,=Y2;U30R20D4L16D22R12U4L4U4R8" 
DRAW "D12L20BE1P=FILL; ,=EDGE;" 

' WTTW 

X2= X2 + 26 

DRAW "C=EDGE ; BM=X2; ,=Y2 ;U30R4D26R12U26R4D30L20BE1P= 
FILL; ,=EDGE;" 

' "RH 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2;,=Y2;U30R20D15L12F15L4H15D15L4BUl 
9BR4U7" 

DRAW "R12D7L12BL1P=FILL; ,=EDGE;" 

f HAN 

X2= X2 + 30 

DRAW "C=EDGE ; BM=X2; ,=Y2;U30R20D30L4U15L12D15L4BU19BR 
auge 

DRAW "R12D7L12BL1LP=FILL; ,=EDGE;" 

' ит" 

Х2= Х2 26 

DRAW "C=EDGE ; BM=X2; ,=Y2;BU26U4R20D4L8D26L4U26L8BE1P= 
FILL: ,=EDGE:" 

' "T" 

X2= X2 + 26 

DRAW "C=EDGE;BM=X2; ,=Y2;U3O0R4D30L4BE1P=FILL;,=EDGE;" 
' HOW 

X2= X2 + 14 

DRAW "C=EDGE ; BM=X2; ,=Y2;U30R20D30L20BR4BU4U22R12D22L 
12BG1P=FILL; ,=EDGE;" 

' UN" 

X2= X2 + 26 

DRAW "C=EDGE ; BM=X2; ,=Y2;U30R4F12U12R4D30L4H12D12L4BE 
LP=FILL; ,=EDGE;" 

'NEW LINE 

t "AH 

DRAW "C=EDGE ;BM=X3 ; ,=Y3;U30R20D30L4U015L12D15L4BU19BR 
4U7" 

DRAV "RI2D7LI2BLIP-FILL, , =EDGE;" 

t "N" 

X3= X3 + 26 

DRAW "C=EDGE;BM=X3; ,=Y3 ;U30R4F12U12R4D30L4H12D12L4BE 
1P=FILL; ,=EDGE;" 

' иди 

X3= X3 + 26 

DRAW "C=EDGE;BM=X3; ,=Y3;U3O0R20D30L4U15L12D15L4BUL9BR 
207" 

DRAW "R12D7L12BL1P=FILL; ,=EDGE;" 

' "r" 

X3= X3 + 26 

DRAW "C=EDGE ; BM=X3; ,=Y3;U30R4D26R16D4L20BE1P=FILL; ,= 
EDGE:" 
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1040 
1950 
1060 
1070 
1080 
1090 
1100 
TETO 
11209 
1130 
1140 
1150 
bY 
DZE 
1180 
1190 
1200 
1210 
1220 
1200 
1240 
1250 
1260 
1270 
1280 
1290 


1300 


1510 


13290 


155 
1340 


1559 


1360 


1979 


1380 


2890 


1400 


1410 


t LAL 

X3= X3 + 26 

DRAW "C=EDGE;BM=X3; ,=Y3;BU30BL1R5F9R4E9R5G14D16L4" 
DRAW "U16H14BR2BD1P=FILL; ,=EDGE;" 

' "su 

X3= X3 + 40 

DRAW "C=EDGE;BM=X3;,=Y3;BUlU2E2R16H19U4E2R20F2D2G2" 
DRAW "L16F19D4G2L20H2BR2BU1P=FILL; , =EDGE;" 

t или 

X3= X3 + 32 

DRAW "C=EDGE;BM=X3;,=Y3;U30R4D30L4BE1P=FILL;,=EDGE;" 
t ис̧и 

X3= X3 + 14 

DRAW "C=EDGE;BM=X3;,=Y3;BUlU2E2R16H19U4E2R20F2D2G2" 
DRAW "L16F19D4G2L20H2BR2BU1P=FILL; , =EDGE;" 

NEXT I 

LOCATE POM 

INPUT "WOULD YOU LIKE INSTRUCTIONS (Y/N)";Z$ 


IF Z$ = "Y" OR ZŚ = "y" THEN GOTO 1240 
GOTO 2000 

SCREEN 0 

VIDTH 80 

COLOR 15,1 

CLS 

PRINT 


PRINT TAB(15) "THIS 
WILL ALLOW YOU TO" 
PRINT TAB(10)"EVALUATE DIFFERENT HELICOPTER CONFIGU 
RATIONS AND DETERMINE" 

PRINT TAB(10)"THEIR LIFT AND DRAG COEFFICIENTS AND 
THEIR EQUIVALENT FLAT" 
PRINT TAB(10) "PLATE AREA. 
THE TUNNEL RUNS WILL" 
PRINT TAB(10) "BE USED WITH THIS PROGRAM." 

PRINT TAB(15) "THE FIRST OPTION IS USED TO RECORD 
THE DATA DURING" 

PRINT TAB(10) "THE TUNNEL RUNS. 
CONVERTS THE RECORDED" 

PRINT TAB(10)"RAW COUNTS TO FORCES AND MOMENTS. 
WILL HAVE THE CHANCE" 

PRINT TAB(10)"TO PROVIDE NAMES FOR EACH OF THESE 
FILES AS YOU GO ALONG." 

PRINT TAB(10)"THE THIRD OPTION LETS YOU EXAMINE THE 
FORCES RECORDED. IT" 

PRINT TAB(10)"WILL BE OF NO USE TO YOU TO EXAMINE THE 
RAW COUNTS." 

PRINT TAB(15)"THE FOURTH OPTION USES THE CONVERTED 
DATA TO COMPUTE" 

PRINT TAB(10)"THE DESIRED PARAMETERS. 
MUST BE CONVERTED USING" 


IS A MENU DRIVEN PROGRAM THAT 


THE DATA GENERATED DURING 


THE SECOND OPTION 


YOU 


THE DATA 
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1420 


1430 


1440 


1450 


1455 
1456 
1460 


1470 


1480 


1490 


1500 


1510 


1520 
"8520 
2000 
2010 
2020 


2030 
2040 
2050 


2070 
2080 
2090 
2100 
2170 
2120 
2030 
2140 
2450 
2160 


2170 
2180 
2100 
2200 


2210 


PRINT TAB(10)"OPTION TWO BEFORE THE PARAMETERS CAN BE 
CALCULATED. YOU WILL" 

PRINT TAB(10)"THEN PROVIDE A NAME FOR THE 
STORE THE CALCULATED" 
PRINT TAB(10) "PARAMETERS. 
THE FIFTH OPTION LETS" 
PRINT TAB(10)"YOU EXAMINE THE DATA THAT CAN BE 
PLOTTED." 
PRINT: INPUT 
CLS 

PRINT TAB(15)"OPTION SIX MUST BE USED TO SORT THE 
X VALUES BEFORE" | 

PRINT TAB(10)"THEY CAN BE PLOTTED. 
LETS YOU ADD A GRID" 

PRINT TAB(10)"TO THE PLOT AND ALSO MARK THE DATA 
Porites. THE PLOTTING IS" 

PRINT TAB(10)"IS EXECUTED WITH THE EIGHTH OPTION. 
YOU CAN PLOT UP TO" 

PRINT TAB(10)"THREE CURVES PER PLOT AND RECEIVE A 
HARD COPY BY PRESSING" 

PRO) "THE CTEL-PRTSC BUTTONS. 
OPTION TERMINATES THIS" 

PRINT TAB(10) "PROGRAM." 

"PRESS RETURN TO CONTINUE";N 


FILES THAT 


USING THESE FILE NAMES, 
"PRESS RETURN TO CONTINUE" ;N 


THE SEVENTH OPTION 


THE FINAL 


SCREEN 0: WIDTH 80: COLOR 14,0,7: CLS: LOCATE 10,1 
PRINT "WHICH DISK DRIVE DO YOU WANT TO STORE YOUR 
DATA FILES ON"; 

PRINT 

PRINT "DRIVE A IS EITHER ON YOUR LEFT OR ON TOP" 
PRINT "DRIVE B IS EITHER ON YOUR RIGHT OR ON THE 
BOTTOM" 

PRINT | | 
INPUT "ENTER [ A OR B ]";FD$: FD$=FD$+":" 
GRD$="NO GRID": MRKS="MARK fl 
OPT$-"CURVE" 

GOSUB 2230 

GOSUB 2400 

X%=5+DX%: Y%=12: LOCATE X%,Y%,1,0,7 
FIRST%=5: LAST%=13 

ANS$=INKEY$: IF ANS$="" THEN 2150 

IF ANS$=CHR$ (0)+CHR$(80) THEN IF X%<LAST% 
X%+1 ELSE X%=FIRST% 

IF ANS$=CHR$(0)+CHR$(72) THEN IF X%>FIRST% THEN X%= 
X%-1 ELSE X%=LAST% 

LOCATE X%,Y%: IF ANS$<>CHR$(13) THEN 2150 
FLAG1%=CSRLIN: FLAG1%=FLAG1%-4 

ON FLAG1% GOSUB 3120,2740,2630,2800,2680,3170,2870, 
3550,3680 

IF FLAG1%<>1 THEN GOTO 2120 ELSE GOTO 2130 


PD$="C:" 


THEN X%= 
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AR と INITIALIZING CONSTANTS---------------- 


2230LM X(300)7/5(300)), C(2700)) 5 (6000 TED 
2240 DIM FILESTK$(10) 
2250 OUT 985,6 
2260 SX%=5: SY%=55 
2270 M$="COMMAND :" 
2280 BLANKS=" i 
2290 B$=" H 
2300 SCREEN 0,1: CLS 
2310 HDPOS%=48 
2320 KEY(1) ON: ON KEY (1) GOSUB 2350 
2330 RETURN 
2340 "------------- CHANGE FOREGROUND COLORS a ae 
2350 FG%=FG%+1 
2360 IF FG%>15 THEN FG%=1 
2370 OUT 985, FG% 
2380 RETURN 
2390 "--------- PRINT MAIN OPTION MENU een --- 
2400 SCREEN 0,1: COLOR 14,0,7: CLS: LOCATE 3,1 
2410 PRINT TAB(5) "xxxxxxxx MAIN OPTIONS MENU 
*kkkkkkkk!! 
2420 PRINT 
2430 PRINT TAB(5) "* _ RECORD TUNNNEL DATA 
JH 
2440 PRINT TAB(5) "* _ REDUCE RAW TUNNEL DATA 
ję 1 
2450 PRINT TAB(5) "* — EXAMINE DATA FILE FOR TUNNEL 
DATA x" 
2460 PRINT TAB(5) "* _ CALCULATE Cd, CI, AND ESF PAR 
x" 
2470 PRINT TAB(5) "x _ EXAMINE DATA FILE FOR PLOTTI 
NG x" 
2480 PRINT TAB(5) "x _ SORT X VALUES 
ki 
2490 PRINT TAB(5) "x _ PLOTTING OPTION MENU 
Јен 
2500 PRINT TAB(5) "x _ DO PLOTTING 
ki 
2510 PRINT TAB(5) "x _ EXIT 
ди 
2520 РЕТМТ 
2530 PRINT TAB(5) ik. s. s. s. kk. 
kkkkkkkk!! 
2540 LOCATE3 5 Era ao. 
2550 LOCATE 4,55 : PRINT "------------------------ f! 
2560 LOCATE SX%,SY% : PRINT "DATA FILE DRIVE = "+FD$ 
2570 LOCATE SX%+2,SY% : PRINT "USER OPTIONS ca 
2580 LOCATE SX%+3,SY% : PRINT OPT$ 
2590 LOCATE SX%+4,SY% : PRINT GRDS$ 
2600 LOCATE SX%+5,SY% : PRINT MRKS 
Figure A.9 MAIN.BAS (cont.) 


40 


2640 
2620 
2630 
2640 
2090 
2660 
2670 
2680 
2550 
2700 
2710 
2720 
203 0 
2740 
2750 
2760 
2 740 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 


2890 
2900 


2910 


2920 


25510 


2940 


2950 
2960 


2270 
2980 
2990 
3000 


3010 


RETURN 

"---DATA EDITOR FOR TUNNEL DATA------------- 
CHAINFILES = PDS+"ADATA. BAS" 

DX%= 

COMMON FD$,PD$,DX% 

CHAIN CHAINFILES 

RETURN 

'=---DATA EDITOR FOR PLOTTING DATA---------- 
CHAINFILES = PD$+"BDATA. BAS" 

DX%= 

COMMON FDS, PDS, Dx% 

CHAIN CHAINFILES 

RETURN 

'----DATA REDUCTION------------------------- 
CHAINFILES = PDS+"RED.BAS" 

DX も =3 

COMMON PDS, FDS, DX る 

CHAIN CHAINFILES 

RETURN 

.----- CALCULATE ROUTINE-------------- 
CHAINFILES = PDS + "COMP. BAS" 

DX%=4 

COMMON PDS, FDS, Dx% 

CHAIN CHAINFILES 


REUTRN 

78” .:c .ETON VEA == 

"1 SCREEN 0,1: LOCATE 6,1 

PRINT TAB(20) "xxxxxx PLOTTING OPTION MENU kkk 

kkk! | 

PRINT 

PRINT TAB(20) "x GRID 6 MARK 
хә 

PRINT TAB(200 "x GRID & NO MARK 
хи 

PRINT TAB(20) "x NO GRID 8 MARK 
kü 

PRINT TAB(20) "x NO GRID & NO MARK 
kn 

PRINT TAB(20) "x EXIT 
je 

PRINT 

PRINT TAB(20) "**kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

kkk 

X%=8: Y%=27: LOCATE X%,Y%,1,0,7 

FIRST%=8: LAST%=12 


ANS$-INKEY$: IF ANSŞ-"" THEN 2990 

IF ANS$=CHRS(0)+CHR$(80) THEN IF X3<LAST% THEN X%= 
X%+1 ELSE X%=FIRST% 

IF ANSS=CHRS (0) +CHRS (72) 
X$+1 ELSE X%=LAST% 


THEN IF X%>FIRSTS THEN X%= 
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3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3145 
3150 
3155 
3160 
21570 


3180 
3190 
3200 
3210 
3220 
3230 
3240 


3250 
3260 
3270 
3280 


3290 
3300 
3310 
3320 
3330 


3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 


LOCATE X%,Y%: IF ANSS<>CHRS(13) THEN 2990 
FLAG%=CSRLIN: FLAG%=FLAG%-7 


ON FLAG$ GOSUB 3070,3080,3090,3100,3110 


DX%=7 

RETURN 

GRDS="GRID 1: MRKS="MARK "RETURN 
GRDS="GRID 1: MRK$="NO MARK" : RETURN 
GRD$="NO GRID": MRK$="MARK " 2 RETURN 


GRD$="NO GRID": 
RETURN 

! ーーーーーー RECORD 
CHAINFILEŞ = 
DX%=1 

COMMON PDS, FDS,DX% 
CHAIN CHAINFILES 
RETURN 


MRKS="NO MARK" RETURN 


TUNNEL DATA--------------- 
PDS + "RUNS. BAS!" 


CLS: INPUT “ENTER NAME OF FILE 
FILES=FDS+FILES 


TO BE SORTED" ELE SE 


PRINT "READING FILE "+FILES+"...": BEEP 
GOSUB 3370 
PRINT "SORTING" 
L%=2: K%=NOD%-1: R%¥=NOD% 
WHILE (L%<=R%) l 
FOR J%=R% TO L% STEP -1 


IF (X(J%-1)>X(J%)) THEN SWAP X(J%),X(J%-1): SWAP 
Y ( け も ) , Y (ける -1 ) : K も = ける 
NEXT 
L%=K%-1 
FOR J% = L% TO R% 
IF (X(J%$-1)>X(J%)) THEN SWAP X(J%) ,X(J%-1): SWAP 
Y(J$),Y(98-1): K$=J% 
NEXT 
R%=K%+1 
WEND 
CLS : BEEP 
PRINT "SORTED FILE ";FILE$;" BEING SAVED ...": GOSUB 
3460 
DX%=6 
RETURN 


OPEN FILES FOR INPUT AS #2 
INPUT #2,NOD%,Y1$,X1$, CONFIGS, CONF 
INPUT 42,0 
FOR J = 1 TO NOD% 
INPUT 42, Y (DEAT) 
NEXT 
CLOSE #2 
RETURN 


Figure A.9 MAIN.BAS (cont.) 


42 


3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 


2260 
DIO 
3580 


3590 
3600 
3601 


3602 


3610 
3620 
5630 
3640 
3650 
3660 
3670 
3680 
3690 


3700 


“sizo 


OPEN FILES FOR OUTPUT AS #1 
WRITE #1,NOD%,Y1$,X1$, CONFIGS, CONF 
WRITE 41,0 | 
FOR J = 1 TO NOD% 
WRITE $1,Y(J),X(J) 
NEXT 
CLOSE #1 
RETURN 
.--------- INITIALIZING PLOTTER----------------- 


CLS: INPUT "NUMBER OF DATA FILE(S) TO BE PLOTTED ON 
THE SAME PLOT =";NOF% 
PRINT "ENTER FILE NAME (S) :":BEEP 
FOR I$ = 1 TO NOF% 
PRINT "FILE 4"; :PRINT USING "t4";1%; INPUT " ="; 
FILEŞ 
FILES = FDS + FILES 
FILESTKS (1%) = FILES 
PRINT "HOW DO YOU WANT THIS CURVE PLOTTED ? ENTER 
(EO SROSJUSTSPLOTU 


INPUT "THE POINTS OR (CURVE) TO PRODUCE A CURVE 
FIT":OPT$ 
CURVES (I%) = 
NEXT 
PFILES = 
SL% = 100 
CHAIN FILE2$,30,ALL 
DX%=8 
RETURN 
CLS: 
LOCATE 10,5: PRINT "PROGRAM TERMINATED, REMOVE DATA 
FILES FROM DISK DRIVE" 


OPTS 


PD$+"MAIN.BAS": FILE2$= PDS + "PLOTTER. BAS" 


LOCATE 11,5: PRINT "AND SECURE COPMUTER AND AMPLIFIE 
GĦ | 
END 
Figure A.9 MAIN.BAS (cont.) 
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190 
200 
277 
220 
230 
240 
250 
260 
270 
280 
290 
300 
210 
. 320 
550 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 


460 
470 
480 


"NAME: 


Data Acquisition And Control (DAAC) 
HEADER for BASICA 


'FILE NAME: DACHDR. BAS 


ee DEVICE NAME: DAAC 


"RESERVED FUNCTION NAMES: 


AINM, AINS, AINSC, AOUM, AOUS, 
BINM, BINS, BITINS, BITOUS, BOUM, BOUS, 
CINM, CINS, CSET, DELAY 


-.. DEF SEG VALUE NAME: DSEG 


'NAMES DEFINED AND USED BY HEADER: 


ADAPT%, AI, COUNT, FOUND% 
HNAMES, SG%, STATS 


'When using the BASICA Interpreter, this header 

'must be executed before any function calls are 

'made that access the DAAC adapter. It initializes 

"a number of variables for each function call. These 
'variables are reserved and should not be used except 
"to access the DAAC adapter. This routine also does a 
'DEF SEG to the segment where the DAAC Device Driver 
'(DAC.COM) is loaded. If you execute a DEF SEG to 
'access other hardware, you must DEF SEG to the segment 
'of the DAAC Device Driver before any subsequent 
‘calls to access the DAAC adapter. 


FOUND% 


SG% 


= 0 
&H2E 


"Start searching the interrupt vectors until you find 
'one that points to the DAAC device driver. 

'Do a DEF SEG to that segment. 

WHILE ((SG% <= &H3E) AND (FOUND% = 0)) 


WEND 


DEF SEG = O 
DSEG = PEEK(SG%) + PEEK(SG% + 1) * 256 
DEF SEG = DSEG 
HNAMES=" " 
FOR AI=10 TO 17 
HNAMES = HNAMES 4 CHRS(PEEK(AI)) 
NEXT AI 
IF HNAMES = "DAAC " AND PEEK(18) + PEEK(19) <> 
O THEN FOUND% = 1 
SG% = SG% + 4 


IF FOUND = O THEN PRINT "ERROR: DEVICE DRIVER DAC.COM 
NOT FOUND" : END 
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490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
780 
740 
750 


760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 


255 


"Nov initialize all function name variables for calls 


'to access the device driver. 

AINM = PEEK(&H13) * 256 + PEEK(&H12) 
AINS = PEEK(&H15) * 256 + PEEK(&H14) 
AINSC = PEEK(&H17) * 256 + PEEK(&H16) 
AOUM = PEEK(&H19) * 256 + PEEK(&H18) 
AOUS = PEEK(&H1B) * 256 + PEEK(&H1A) 
BINM = PEEK(&H1D) * 256 + PEEK(&H1C) 
BINS = PEEK(&H1F) * 256 + PEEK(&H1E) 
BITINS = PEEK(&H21) * 256 + PEEK(&H20) 
BITOUS = PEEK(&H23) * 256 + PEEK(&H22) 
BOUM = PEEK(&H25) * 256 + PEEK(&H24) 
BOUS = PEEK(&H27) * 256 + PEEK(&H26) 
CINM = PEEK(&H29) * 256 + PEEK(&H28) 
CINS = PEEK(&H2B) * 256 + PEEK(&H2A) 
CSET = PEEK(&H2D) * 256 + PEEK(&H2C) 
DELAY = PEEK(&H2F) * 256 + PEEK(&H2E) 
'Finally, execute any call to re-initialize the 
'device driver from any former invocation of BASIC. 
ADAPT% = 0 

COUNT = 1 

STATS = 0 

CALL DELAY (ADAPT%, COUNT, STATS) 


'End of DAAC BASICA Header 
' 
КЕМ-- RUNS.BAS; 
WIND 
REM-- TUNNEL RUNS 
KEY FF: COLOR 25,1,4: 
GES 
PRINT: PRINT 
PRINT TAB(10) 


PROGRAM TO RECORD THE DATA FROM 


CLS 


"HELO CONFIGURATIONS" 


PRIN@: PRINT 

PRINT TAB(5) "1. ATTACK NOSE, STRAIGHT TAIL" 
PRINT TAB(5) "2. ATTACK NOSE, LOW TAIL" 
PRINT TAB(5) "3. ATTACK NOSE, HIGH TAIL" 
PRINT TAB(5) "4. SMOOTH NOSE, STRAIGHT TAIL" 
PRIN TAB(59 "5. SMOOTH NOSE, LOW TAIL" 
ERINT TAB(5) "6. ISMOOTH NOSE, HIGH TAIL" 
EELUU TAB(5) "7. BLUNT N@SE, STRAIGHT TAIL" 
PRINT TAB(5) "8. BLUNT NOSE, LOW TAIL" 
PRAMSTAB(5) "9. BLUNT NOSE, HIGH TAIL" 
PRINT: PRINT 

INPUT "WHICH CONFIGURATION IS BEING RUN";N 
PRINT 

INPUT "DOES CONFIGURATION INCLUDE LANDING GEAR";Z$ 


ON N GOSUB 4550,4620,4690,4750,4820,4870,4920,4990, 


5040 


Figure A.10 
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WE TGT TS ARE T INI PORNDS=--=————- 


RUNS.BAS (cont.) 


960 

970 

980 

990 

1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
14750 
1130 
1140 
1150 
1160 


1170 
1180 
1190 
1200 
1210 
2670 


2671 


2672 


2673 
2690 


2700 
2710 
2720 


2730 
2740 


2750 
2760 
27: 
2780 
2790 
2800 
2810 


IF CONF = 10 THEN W = 22.7892 
IF CONF = 15 THEN W = 23.4623 
IF CONF = 20 THEN W = 24.6945 
IF CONF = 25) THEN W = 25.2670 
IF CONF = 30 THEN W = 24.6945 
IF CONE =35) THEN W = 2552676 
IF CONF = 40 THEN W = 20.8134 
IF CONF = 45 THEN W = 21.5866 
IF CONF =150 TAEN W = 22:7137 
IF CONF = 55 THEN W = 23.4913 
IF CONF = 60 THEN W = 22.7187 
IF CONF = 65 THEN W = 23.4913 
IF CONF = 70 THEN W = 21.3315 
IF CONF = 75 THEN W = 23.2892 
IF CONF = 80 THEN W = 23.2368 
IF CONF = 85 THEN W = 25.1945 
IF CONF = 90 THEN W = 23.2368 
IF CONF = 95 THEN W = 25.1945 
CLS | 
КВИ-= ә RECORD TUNNEL- DATA- ни 


DIM L(100) ,D(100) ,Y(100) ,PM(100) ,YM(100) ,RM(100) ,AOA( 
100) 
DIM DAT(399) ,DAT%(399) ,DAT1(399) ,DAT1%(399) 


COLOR 15,1::5570EE: 78055 

PRINT "RECORDING OF VIND TUNNEL RAV DATA" 

INPUT "WHAT IS THE TUNNEL SPEED (Q) FOR THIS RUN";Q > 
PRINT 

PRINT 'WITHOUT THE MODEL IN PLACE, ADJUST THE ZERO 


SET SCREWS TO ZERO" 

PRINT “OUT EACH AMPLIFIER. 
ZEROED, LOAD THE" 

PRINT "MODEL ONTO THE STING. 
MOUNTED, PRESS RETURN" 


AFTER ALL AMPLIFIERS ARE 


WHEN THE MODEL IS 


INPUT "TO RECORD A NO FORCE ZERO READING";N 

PRINT: PRINT "n ZD ZL ZY ZPM ZYM 
ZRM ZAOA" 

STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 0: DEVICES= 9: CHANLO%= O: CTRL$= 0: CHANHI%= 


CALL AINSC(ADAPT%, DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 
,STORS, COUNT, RATE, DAT% (0) , STATS) 

20 = 0:2РМ=0:21=0:2УМ=0 

IF STAT%<> О ТНЕМ РЕТМТ USING "EXECUTION ERROR ###"; 
STATS: END 

FOR J = 0 TO 396 STEP 4 

DAT (J) =(DAT%(J)/204.8)-10 

ZD =ZD + DAT(J) 


NEXT J 
ZD = ZD/100 
FOR J = 1 TO 397 STEP 4 


DAT(J)=(DAT%(J)/204.8)-10 
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2820 
2830 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
20 
2930 
2940 


2950 
2960 


2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3051 
3052 
3053 
3054 
3060 


3070 


3080 
3090 
2100 
3110 
3 MEL 
3120 
3530 
3140 


3450 
3160 


3170 
3 130 
5 EDU 


ZL = ZL + DAT(J) 


NEXT J 
FOR J = 2 TO 398 STEP 4 
DAT (J) =(DAT% (J) /204.8)-10 


ZPM = ZPM + DAT(J) 

NEXT J 

FOR J = 3 TO 399 STEP 4 

DAT(J)=(DAT%(J)/204.8)-10 

ZYM = ZYM + DAT(J) 

NEXT J 

STATS=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPTS= 1: DEVICES= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=2 
CALL AINSC (ADAPT$, DEVICE% , CHANLO% , CHANHI% , CTRL$ , MODE% 


, STOR% , COUNT, RATE,DAT1% (0) , STATS) 
ZY =0:ZRM=0:ZAOA=0 

IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

FƏR J = 0 TO 297 STEP 3 

DAT1 (J) =(DAT1% (J) /204.8) -10 

ZRM = ZRM + DAT1(J) 

NEXT 

ZRM = ZRM/100 

FOR J = 1 TO 298 STEP 3 

DAT1 (J) =(DAT1% (J) /204.8)-10 

w = ZY ОАТ1(Ј) 

NEXT J 


FOR J = 2 TO 299 STEP 3 

DAT1(J)=(DAT1%(J)/204.8)-10 

ZAOA = ZAOA + DAT1(J) 

NEXT J 

ZL=ZL/100: ZPM=ZPM/100:ZYM=ZYM/100:ZY=ZY/100:ZAOA=ZAOA 


/ 160 


WOCAPĘ 11,1: PRINT USING "+4. 47#";ZD: LOCATE 11,10: 
PRINT USING "+#.###";7ZL 

LOCATE 11,19: PRINT USING "+#.###";ZY 

LOCATE 11,28: PRINT USING "+#.###";ZPM 

LOCATE 11,37: PRINT USING "+#.###";ZYM 

LOCATE 11,46: PRINT USING "+#.###"; ZRM 

LOCATE 11,55: PRINT USING "+#.###";ZAOA 

REM 

PRINT: PRINT 

INPUT "AFTER PLACING ALL CAL SWITCHES TO + SETTING HIT 
RETURN"; X 

PRINT 

PRINT " CALD CLL CALY CALPM CALYM 
CALRM CALAOA" 

STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 
CALL AINSC(ADAPT%, DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 


, STOR る , COUNT, RATE, DAT% (0) , STATS) 


Figure A.10 RUNS.BAS (cont.) 
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3200 
3210 


3220 
3230 
3240 
3250 
3260 
3270 
3230 
3290 
3300 
3 JO 
3320 
3330 
3340 
3550 
3360 
3370 
3380 
3390 
3400 
3410 


3420 
3430 


3440 
3450 
3460 
3470 
3480 
3490 
3500 
35190 
3520 
3521 
7 
3525 
3524 
2530 


3540 
3550 
3560 
3570 
3580 
3590 
2501 
5592 


CALD=0 : CLL=0 : CALYM=0 : CALPM=0 

IF STAT%<> 0,THEN PRINT USING "EXECUTION ERROR ###" 
; STAT% : END 

FOR J = 0 TO 396 STEP 4 

DAT (J) =(DAT% (J) /204.8)-10 

CALD =CALD + DAT(J) 


NEXT J 
CALD = CALD/100 
1 TO 397 STEP 4 


FOR J = 
DAT (J) =(DAT% (J) /204.8)-10 


CLL = CLL + DAT(J) 
NEXT J 
FOR J = 2 TO 398 STEP 4 


DAT (J) =(DAT% (J) /204.8)-10 


CALPM = CALPM + DAT(J) 
NEXT J 
FOR J = 3 TO 399 STEP 4 


DAT(J)=(DAT$%(J)/204.8)-10 


CALYM = CALYM + DAT(J) 

NEXT J 

STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 1: DEVICE%= CHANLO%= 0: CTRL%= 0: CHANHI%=2 


CALL AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL$%,MODE% 
, STOR% , COUNT, RATE, DATI% (0) STATS) 

CALY ZO CA ERM= OGA AA Ы 

IF STAT%<> O THEN PRINT USING "EXECUTION ERROR ###"; 
STAT%: END 

FOR J = O TO 297 STEP 3 

DAT1(J)=(DAT1%(J)/204.8)-10 


CALRM = CALRM + DAT1(J) 
NEXT J 

CALRM = CALRM/100 

FOR J = 1 TO 298 STEP 3 


DAT1(J)=(DAT1%(J)/204.8)-10 

CALY = CALY + DAT1(J) 

NEXT J 

FOR J = 2 TO 299 STEP 3 

DAT1(J)=(DAT1%(J)/204.8)-10 

CALAOA = CALAOA + DAT1(J) 

NEXT J 

CLL=CLL/ 100: CALPM=CALPM/ 100: CALYM=CALYM/100: CALY=CALY 
/100:CALAOA=CALAOA/100 


LOCATE 18,1: PRINT USING "+#.###";CALD 
LOCATE 18,10: PRINT USING "+#.###";CLL 
LOCATE 18,19: PRINT USING "+#.###";CALY 
LOCATE 18,28: PRINT USING "+#.###" ;CALPM 
LOCATE 18,37: PRINT USING "+#.###";CALYM 
LOCATE 18,46: PRINT USING "+#.###"; CALRM 
LOCATE 18,55: PRINT USING "+#.###";CALAOA 
INPUT "PRESS" RETURN LTOSCOHTELUE” 1 


Figure A.10 RUNS.BAS (cont.) 
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3593 
3600 


3610 
3620 
3630 
3640 


3650 
5651 
3652 


3653 
3654 
3659 
3656 


2037 
3660 
3670 
3680 


3690 


2700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
2510 
3820 
3830 
3840 
3850 
3860 
3870 
3880 


CLS 

PRINT: PRINT "REPLACE THE CAL SWITCHES TO THE CENTER 
POSITION" 

PRINT "AFTER YOU HAVE GOTTEN THE WIND TUNNEL UP TO 


SPEED AND ARE READY" 
PRINT "TO RECORD DATA PRESS RETURN. 
ARE PRINTED THE" 


ONCE THE HEADINGS 


PRINT "F2 KEY WILL RECORD THE DATA. THE Fl KEY WILL 
SAVE THE DATA AND " 

PRINT "RERUN YOU TO THE MAIN MENU. YOU CAN TAKE READ 
INGS FOR ANY ANGLE" 

PRINT "OF ATTACK BETWEEN +10 AND -8 DEGREES." 

PRINT 

PRINT "THIS PROGRAM CORRECTS FOR DRIFT IN THE BALANCE 


AND AMPLIFIERS." 

PRINT "YOUR FIRST AND LAST DATA POINTS SHOULD BE TAKEN 
AT ZERO ANGLE OF" 

PRINT "ATTACK. ALSO, DO NOT TAKE MORE THAN ONE DATA 
POINT FOR THE SAME " 


PRINT “ANGLE OF ATTACK EXCEPT FOR YOUR FIRST AND LAST 
DATA POINT. you" 
PRINT "ARE LIMITED TO 100.DATA POINTS. REMEMBER, EACH 
RUN IS FOR ONE" 
PRINT "SPECIFIC 'Q' SETTING" 
INPUT "PRESS RETURN TO CONTINUE";N 
CLS 
PRINT " DRAG LIFT YAW PITCH YAW 
ROLL AOA TUNNEL" 
PRINT " MOM. MOM. 
MOM. SPEED" 
SOAP = 0: N=4 
FOR K = 1 TO 100 
ON KEY(1) GOSUB 3820 "SET STOP FLAG 
ON KEY(2) GOSUB 3850 "RECORD DATA 
KEY(1) ON: KEY(2) ON 
IF SOAP = 2 THEN GOTO 3780 
IF SOAP = 1 THEN GOTO 4370 
GOTO 3720 
SOAP = 0 
NOD% = K 
NEXT K 
GOTO 4370 
REM SET STOP FLAG 
SOAP = 1 
RETURN 
REM STEPS TO RECORD DATA 
STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 0: DEVTCE も = CHANLO%= 0: CTRL%= 0: CHANHT も =3 
CALL AINSC(ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 


, STORS , COUNT, RATE, DAT%(0) , STATS) 


Figure A.10 RUNS.BAS (cont.) 
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3890 
3900 


3910 
3020 
3930 
3940 
3950 
3960 
3970 
3980 
20 
4000 
4010 
4020 
4030 
4040 
4041 
4050 


D(K)=0: L(K) =0: YM(K) =O: PM(K) =0 
IF STAT%<> 0,THEN PRINT USING "EXECUTION ERROR +++"; 


STAT%: END 

FOR J = 0 TO 396 STEP 4 
DAT (J) =(DAT% (J) /204.8)-10 
D(K) =D(K) + DAT(J) 

NEXT J 
D(K) = 
FOR J = 
DAT (J) = 
L(K) = 
NEXT J 
FOR J = 2 TO 398 STEP 4 
DAT (J) =(DAT% (J) /204.8)-10 
PM(K) = PM(K) + DAT(J) 
NEXT J 

FOR J = 3 TO 399 STEP 4 
DAT (J) =(DAT% (J) /204.8)-10 
YM(K) = YM(K) + DAT(J) 


D(K) /100 

1 TO 397 STEP 4 
(DAT% (J) /204.8)-10 
L(K) + DAT(J) 


4060 
4070 
4080 
4090 


NEXT J 
STATS=0: MODE%=0: STOR%=0: COUNT=1007 RATE—500 
ADAPT%= 1: DEVICE%= 9: CHANLOS= 0: CTRLSE 0: CHANHI%=2 
CALL AINSC (ADAPTS, DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 
, STOR , COUNT, RATE, DATI: (0) , STATS) 
Y(K) =0:RM(K)=0: AOA(K)=0 
IF STAT%<> O THEN PRINT USING "EXECUTION ERROR +++"; 
STAT% : END 
FOR J =O TO 297 STER 3 
DAT1(J)=(DAT1%(J)/204.8)-10 
RM(K) = RM(K) + DAT1(J) 
NEXT J 
RM(K) = RM(K)/100 
FOR J = 1 TO 298 STEP 3 
DAT1 (J) =(DAT1%(J)/204.8)-10 
Y(K) = v(K) + DATI(J) 
NEXT J 
FOR J = 2 TO 299 STEP 3 
DAT1/-J)=(DAT1%/(J)/204.8 GS 
AOA(K) = AOA(K) + DAT1(J) 
NEXT J 
AOA(K)=AOA(K)/100:L(K)=L(K)/100:PM(K)=PM(K)/100:YM(K) 
=YM(K)/100:Y(K)=Y(K)/100 
4255 В = АОА(К) – 2АОА 
4256 AOA(K) = (5.8469*B) + (.0077583*(B"2)) 
4260 IF N> 23 THEN N=4:CLS:PRINT " DRAG 
PITCH YAW ROLL AOA 

MOM. MOM. MOM. 
LOCATE N,1l: PRINT USING "+#.###";D(K) 
LOCATE N,10: PRINT USING "+#.###";L(K) 
LOCATE N,19: PRINT USING "+#.###";Y(K) 


4100 
4110 


4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 


CIRE YAW 
TUNNIS SS BE 


SPEED" 
4270 


4280 
4290 
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4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4371 
4372 


2177: 


4374 
4375 


4376 
4377 
4378 
4379 
4380 
4381 
4382 
4383 
4385 
4386 
4387 
4390 
4400 


4410 
4420 
4421 
4430 
4470 
4480 
4490 
4500 
4510 
4520 
4521 


4522 


4523 
4524 
4525 
4530 
4531 
4532 


LOCATE N,28: 
LOCATE N,37: 
LOCATE N,46: 
LOCATE N,55: 
N=N+1 
SOAP = 2 
RETURN 
SOAP = 0 
"---CORRECT FOR DRIFT------------- 
DIFFL =(L(NOD%)-L(1)): DIFFD=(D(NOD%)-D(1)):DIFFY=(Y( 
NOD%) -Y(1)) 
DIFFPM=(PM(NOD%)-PM(1)) :DIFFYM=(YM(NOD%)-YM(1)) : DIFFR 
M= (RM(NOD%) =RM(1) ) : DIFFAOA= (AOA (NOD%) =AOA (1) ) 
D=NOD%-1 
CORL=DIFFL/D:CORD=DIFFD/D:CORY=DIFFY/D:CORPM=DIFFPM/D 
: CORYM=DIFFYM/D:CORRM=DIFFRM/D:CORAOA=DIFFAOA/D 
A=1 
FOR K = 2 TO NOD% 

L(K)=L(K) - (A*CORL) 

D(K)=D(K) - (A*CORD) 

Y(K)=Y(K) - (A*CORY) 

PM (K) =PM(K) - (A*CORPM) 

YM (K) =YM(K) -(A*CORYM) 

RM (K) =RM (K) - (A*CORRM) 

A = A+1 
NEXT K 


PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 


"44,224" 2 PM(K) 
ER. FRE ;YM(K) 
"+4. 441"; RM(K) 
Wee tee" s AOA (K) 


: PRINT 
INPUT "WHAT IS THE NAME FOR THE FILE TO STORE THE RAW 
DATA"; FILES 
FILES-FDS-FILES 
OPEN FILES FOR OUTPUT AS #1 
WRITE #1,NOD%, CONFIGS, CONF 
WRITE #1,9,W 
WRITE +1,ZL,ZD,ZY,ZPM,ZYM, ZRM, ZAOA 
WRITE #1,CLL,CALD, CALY, CALPM, CALYM, CALRM, CALAOA 
FOR J=1 TO NOD% 


WRITE $1,L(J) ,D(J),Y(J),PM(J) ,YM(J) ,RM(J) , AOA (J) 
NEXT J 

CLOSE #1 

CLS:PRINT "DO YOU WANT TO TAKE ANOTHER SET OF DATA 


POINTS FOR THE SAME" 

INPUT "CONFIGURATION BUT A DIFFERENT 'Q' 

) "¡ANSS 

IF ANS$="Y" OR ANS$="y" THEN GOTO 4524 ELSE GOTO 4530 
PRINT: INPUT “WHAT IS THE “Q” VALUBeFOR THIS RUN" ;Q 
GOTO 3610 

!————-RETURN TO MAIN PROGRAM 

COMMON PD$,FD$,DX% 

CHAIN PD$+"MAIN.BAS",2110 


SETTING (Y/N 


Figure A.10 RUNS.BAS (cont.) 
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4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 


REM----SUBROUTINES FOR CONFIGURATION MARKING------- 
IF Z$ = "Y" OR ZS = "y" THEN GOTO 4590 
CONF = 10 


CONFIGS="ATTACK NOSE/STRAIGHT TAIL WITHOUT GEAR" 
GOTO 4610 

CONF = 15 

CINFIGS="ATTACK NOSE/STRAIGHT TAIL WITH GEAR" 
RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4660 

CONF = 20 

CONFIGS="ATTACK NOSE/LOW TAIL WITHOUT GEAR" 

GOTO 4680 l 

CONF € 25 

CINFIGS="ATTACK NOSE/LOW TAIL WITH GEAR" 

RETURN 

IF Z$ = "Y" OR ZS = "y" THEN GOTO 4730 

CONF = 30 

CONFIG$="ATTACK NOSE/HIGH TAIL WITHOUT GEAR" 

GOTO 4740 

CONF = 35: CINFIGS$="ATTACK NOSE/HIGH TAIL WITH GEAR" 
RETURN 

IF Z$ = "Y" OR Z$ = "y" THEN TCOTOTA790 

CONF = 40 

CONFIG$="SMOOTH NOSE/STRAIGHT TAIL WITHOUT GEAR" 
GOTO 4810 

CONF = 45 

CINFIGS="SMOOTH NOSE/STRAIGHT TAIL WITH GEAR" 
RETURN 

IF Z$ = "Y" OR ZS = "y" THEN GOTO 4850 

CONF=50: CONFIGS="SMOOTH NOSE/LOW TAIL WITHOUT GEAR" 
GOTO 4860 

CONF = 55: CINFIGS=""SMOOTH NOSE/LOW TAIL WITH GEAR" 
RETURN f . 

IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4900 
CONF=60:CONFIG$="SMOOTH NOSE/HIGH TAIL WITHOUT GEAR" 
GOTO 4910 

CONF = 65: CINFIG$="SMOOTH NOSE/HIGH TAIL WITH GEAR" 
RETURN | | 
IF Z$ = "Y" OR Z$ = "y" THEN GOTO 4960 

CONF = 70 

CONFIGS="ATTACK NOSE/STRAIGHT TAIL WITHOUT GEAR" 
GOTO 4980 

CONF = 75 

CINFIGS="ATTACK NOSE/STRAIGHT TAIL WITH GEAR" 
RETURN 

ТЕ 25 = "У" OR Z$ = "y" THEN GOTO 5020 

CONF = 80:CONFIGS="BLUNT NOSE/LOW TAIL WITHOUT GEAR" 
GOTO 5030 

CONF = 85: CINFIGS="BLUNT NOSE/LOW TAIL WITH GEAR" 
RETURN 


Figure A.10 RUNS.BAS (cont.) 
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5040 TZS = SOR ZS = "vy" THEN GOTO 5070 

5050 CONF=90:CONFIG$="BLUNT NOSE/HIGH TAIL WITHOUT GEAR" 
5060 GOTO 5080 

5070 CONF = 95: CINFIGS="SMOOT NOSE/HIGH TAIL WITH GEAR" 
5080 RETURN 


Figure A. 10 SRUNS=eAS= (cont. ) 
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217 


250 


260 
270 


280 


290 
300 


310 
320 


330 
340 


350 
360 


REM----RED.BAS (CONVERT RAW COUNTS TO FORCES) ---- 
REM 
REM 
COLOR 15,1:KEY OFF: CLS 
DIM L(100), LI (1009, 12 (1669 NN FONT 
(100) Y111906077 77100) 
DIM PM(100) ,PM1(100) ,PM2(100) ,YM(100) ,YM1(100) ,YM2(10 
0),RM(100) ,RM1(100) ,RM2 (100) 
DIM LF(100) ,DF (100) , AOA (100) 
PRINT: PRINT"WHAT IS THE NAME OF THE FILE THAT CONTAINS 
THE RAW DATA" 
INPUT"THAT YOU WISH TO CONVERT"”FILEŞ 
FILES=FDS+FILES 
.----- READ IN RAV DATA FILE----- 
OPEN FILES FOR INPUT AS fl 
INPUT #1,NOD%, CONFIGS, CONF 
INPUT #1,0,W 
INPUT $Ż1,ZL,ZD,ZY,ZPM,ZYM,ZRM,ZAOA 
INPUT f1,CLL,CALD,CALY,CALPM,CALYM,CALRM,CALAOA 
FOR J = 1 TO NOD% 
INPUT $Ż1,L(J) ,D(J) ,Y(J) ,PM(J) , YM(J) ,RM(J) ,AOA(J) 
NEXT J 
CLOSE fl 
PRINT:PRINT"PLEASE WAIT WHILE PERFORMING CONVERSION" 
PRINT "THE CONVERSION WILL TAKE APPROXIMATELY TEN MIN 
UTES FOR " 
PRINT "TWENTY DATA POINTS" 
.------ READ IN CALIBRATION CONSTANTS-------- 
OPEN "C:CONST" FOR INPUT AS #1 
INPUT #1, INCALL, INCALD, INCALY , INCALPM, INCALYM, INCALRM 
INPUT $41,K1LPOS ,K2LPOS,K1DPOS,K2DPOS,K1YPOS,K2YPOS 
INPUT #1,K1PMPOS , K2PMPOS, K1YMPOS , K2YMPOS ,K1RMPOS, 
K2RMPOS 
INPUT #1,K1LNEG,K2LNEG, K1DNEG, K2DNEG, K1YNEG, K2YNEG 
INPUT 41,К1РММЕС, К2РММЕС, КЈУММЕС,К2УММЕС, КЈЕММЕС, 
К2ЕММЕС 
INPUT 41,DDDLIP,DDDL2P,DVDLIP, DVDL2P,DPMDLIP,DPMDL2P, 
DYMDLIP,DYMDL2P 
INPUT #1,DRMDL1P, DRMDL2P 
INPUT #1,DLDD1P,DLDD2P,DYDD1P, DYDD2P,DPMDD1P,DPMDD2P, 
DYMDD1P, DYMDD2P 
INPUT #1,DRMDD1P, DRMDD2P 
INPUT +41,DLDY1P,DLDY2P,DDDY1P,DDDY2P,DPMDY1P,DPMDY2P, 
DYMDY1P,DYMDY2P 
INPUT #1,DRMDY1P,DRMDY2P 
INPUT #1,DLDPM1P, DLDPM2P,DDDPM1P,DDDPM2P,DYDPM1P, DYDP 
M2P, DYMDPM1P, DYMDPM2P 
INPUT #1,DRMDPM1P, DRMDPM2P 
TNPUT #1 , DLDYM1P , DLDYM2P , DDDYM1P , DDDYM2P , DYDYM1P , DYDY 
M2P , DPMDYM1P , DPMDYM2p 


Figure A.11 RED.BAS - Data Reduction Program 
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370 
380 


390 
400 


410 
420 


430 
440 


450 
460 


470 
480 


490 
500 


510 
520 
590 
540 
550 
560 
570 
580 
590 
600 
610 


Gi っ 


620 


625 
630 


655 
640 


645 
650 


655 
660 


665 
670 


INPUT 41, DRMDYM1P, DRMDYM2P 

INPUT 41,DLDRM1P, DLDRM2P, DDDRM1P, DDDRM2P, DYDRM1P, DYDR 
M2P, DPMDRM1P, DPMDRM2P 

INPUT 41,DYMDRM1P, DYMDRM2P 

INPUT #1,DDDL1N,DDDL2N, DYDL1N, DYDL2N, DPMDL1N, DPMDL2N, 
DYMDL1N , DYMDL2N 

INPUT #1,DRMDL1N, DRMDL2N 

INPUT #1,DLDD1N , DLDD2N , DYDD1N , DYDD2N , DPMDD1N , DPMDD2N , 
DYMDDIN, DYMDD2N | 

INPUT #1,DRMDD1N,DRMDD2N 

INPUT #1,DLDY1N, DLDY2N, DDDY1N, DDDY2N, DPMDY1N,DPMDY2N, 
DYMDY1N, DYMDY2N 

INPUT #1,DRMDY1N,DRMDY2N 

INPUT #1,DLDPM1N, DLDPM2N,DDDPM1N, DDDPM2N, DYDPM1N, DYDP 
M2N, DYMDPM1N 

INPUT £1,DYMDPM2N, DRMDPM1N , DRMDPM2N 

INPUT 41,DLDYM1N, DLDYM2N, DDDYM1N, DDDYM2N, DYDYM1N,DYDY 
M2N,DPMDYMIN 

INPUT #1,DPMDYM2N,DRMDYM1N, DRMDYM2N 

INPUT #1,DLDRM1N,DLDRM2N, DDDRM1N, DDDRM2N, DYDRM1N,DYDR 
M2N,DPMDRMIN 

INPUT #1,DPMDRM2N, DYMDRM1N , DYMDRM2N 

CLOSE #1 

リー ニーーー CONVERT RAW COUNTS TO FORCES------ 

FOR J = 1 TO NOD% 

A = ((INCALL/(CLL-ZL))*(L(J)-ZL)) 


B = ((INCALD/(CALD-2D))*(D(J)-2ZD)) 
C = ((INCALY/ (CALY-ZY) )*(Y¥ (J) -ZY) ) 
О=( (ТМСАГРМ/(САГРМ-2РМ))*(РМ(Ј)-2РМ)) 
E=( (INCALYM/ (CALYM-ZYM) ) * (YM(J) -ZYM) ) 
F=( (INCALRM/ ((CALRM-ZRM) ) * (RM (J) =ZRM) ) 

IF A <= O THEN Kl=KlLNEG:K2=K2LNEG ELSE K1=K1LPOS 
: K2=K2LPOS | 

L(J) = (K1*A) 4 (K2”A”2) 

IF B <= 0 THEN K1=K1DNEG:K2=K2DNEG ELSE K1=K1DPOS 
: K2=K2DPOS 

D(J) = (K1*B) + (K2*B^2) 

IF C <= O THEN K1=K1YNEG:K2=K2YNEG ELSE Kl=K1YPOS 
: K2=K2YPOS 

Y(J) = (Kl*C) + (K2*C^2) 

IF D <= 0 THEN Kl=KLPMNEG:K2=K2PMNEG ELSE Kl=Kl 


PMPOS : K2=K2 PMPOS 

PM(J) = (K1*D) + (K2*D"2) 

IF E <= 0 THEN K1=K1YMNEG:K2=K2YMNEG ELSE K1=Kl 
YMPOS : K2=K2YMPOS 

YM(J) = (K1*E) + (K2*E‘2) 

IF F <= 0 THEN Kl=K1RMNEG:K2=K2RMNEG ELSE Kl=K1 
RMPOS : K2=K2RMPOS 

RM(J) = (K1*F) + (K2*F2) 
NEXT J 
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680 
690 
695 
700 
710 


715 
720 


730 
740 


745 
750 


760 
770 


TAS 
780 


790 
800 


805. 


810 


820 
830 


885 
840 


850 
860 


865 
870 


ニー ニー デー PERFORM INTERACTON CORRECTIONS---- 
FOR K = 1 TO NOD% 


IF L(K) > O THEN GOTO 710 ELSE GOTO 720 

DD1=DDDL1P: DD6=DDDL2P: DY1=DYDL1P: DY6=DYDL2P: DPM1=DP 
MDL1P: DPM6=DPMDL2P: DYM1=DYMDL1P: DYM6=DYMDL2P: DRM1=D 
RMDL1P: DRM6=DRMDL2P 

GOTO 730 
DD1=DDDL1N:DD6=DDDL2N:DY1=DYDL1N:DY6=DYDL2N:DPM1=DP 
MDL1N : DPM6=DPMDL2N : DYM1=DYMDLI1N : DYM6=DYMDL2N : DRM1=D 
RMDLIN: DRM6=DRMDL2N 

IF D(K) > 0 THEN GOTO 740 ELSE GOTO 750 
DL1=DLDD1P:DL6=DLDD2P:DY2=DYDD1P:DY7=DYDD2P:DPM2=DP 
MDD1P: DPM7=DPMDD2 P: DYM2=DYMDD1P: DYM7=DYMDD2P: DRM2=D 
RMDD1P: DRM7=DRMDD2P 

GOTO 760 

DL1=DLDD1N: DL6=DLDD2N: DY2=DYDD1N: DY7=DYDD2N : DPM2=DP 
MDDIN : DPM7=DPMDD2N : DYM2=DYMDDIN: DYM7=DYMDD2N : DRM2=D 
RMDD1N: DRM7=DRMDD2N 

IF Y(K) > 0 THEN GOTO 770 ELSE GOTO 780 

DL2=DLDY1P: DL7=DLDY2P: DD2=DDDY1P: DD7=DDDY2P: DPM3=DP 
MDY1P:DPM8=DPMDY2P:DYM3=DYMDY1P:DYM8=DYMDY2P:DRM3=D 
RMDY LP: DRM8=DRMDY2P 

GOTO 790 

DL2=DLDY1N: DL7=DLDY2N: DD2=DDDY1N: DD7=DDDY2N: DPM3=DP 
MDY1N:DPM8=DPMDY2N:DYM3=DYMDY1N:DYM8=DYMDY2N:DRM3=D 
RMDY1N: DRM8=DRMDY2N 

IF PM(K) > 0 THEN GOTO 800 ELSE GOTO 810 
DL3=DLDPM1P: DL8=DLDPM1P: DD3=DDDPM1P: DD8=DDDPM2P: DY3 
=DYDPM1P: DY8=DYDPM2P: DYM4=DYMDPM1P: DYM9=DYMDPM2P:DR 
M4=DRMDPM1P: DRM9=DRMDPM2P 

GOTO 820 

DL3=DLDPM1N: DL8=DLDPM1N: DD3=DDDPM1N: DD8=DDDPM2N: DY3 
=DYDPM1N: DY8=DYDPM2N : DYM4=DYMDPMIN : DYM9=DYMDPM2N : DR 
M4=DRMDPM1N: DRM9=DRMDPM2N 

IF YM(K) > O THEN GOTO 830 ELSE GOTO 840 
DL4=DLDYM1P:DL9=DLDYM2P:DD4=DDDYM1P:DD9=DDDYM2P:DY4 
=DYDYM1P: DY9=DYDYM2P: DPM4=DPMDYM1P: DPM9=DPMDYM2P:DR 
M5=DRMDYM1P : DRM10=DRMDYM2P 

GOTO 850 
DL4=DLDYM1N:DL9=DLDYM2N:DD4=DDDYM1N:DD9=DDDYM2N:DY4 
=DYDYM1N : DY9=DYDYM2N: DPM4=DPMDYM1N : DPM9=DPMDYM2N: DR 
M5=DRMDYM1N : DRM1O=DRMDYM2N 

IF RM(K) > 0 THEN GOTO 860 ELSE GOTO 876 
DL5=DLDRM1P:DL10=DLDRM2P:DD5=DDDRM1P:DD10=DDDRM2P:D 
Y5=DYDRM1P:DY10=DYDRM2P: DPM5=DPMDRM1P : DPM10=DPMDRM2 
P:DYM5=DYMDRM1P:DYM10=DYMDRM2P 

GOTO 880 

DL5=DLDRMIN : DL10=DLDRM2N : DD5=DDDRMIN : DD10=DDDRM2N:D 
Y5=DYDRM1N : DYLO=DYDRM2N : DPM5=DPMDRMIN : DPM10=DPMDRM2 
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890 
900 
950 
960 
970 


215 
980 


990 
1000 


1005 
TORO 


1020 
1030 


1955 
1040 


1050 
1060 


1070 
1080 


1090 
1100 


1405 
4410 


1120 
1945 0 


15 5 


N: DYM5=DYMDRM1N : DYM10=DYMDRM2N 


L1 (K)=L(K)-(DL1*D(K))-(DL6*(D(K)^2)) 
D1 (K)=D(K)-(DD1*L(K))-(DD6*(D(K)^2)) 


----- SECOND CONSTANT DETERMINATION----------- 


IF L1(K) > O THEN GOTO 970 ELSE GOTO 980 
DD1=DDDL1P:DD6=DDDL2P:DY1=DYDL1P:DY6=DYDL2P:DPM1=DP 
MDL1P: DPM6=DPMDL2P: DYM1=DYMDL1P: DYM6=DYMDL2P: DRM1=D 
RMDL1P: DRM6=DRMDL2P 

GOTO 990 
DD1=DDDL1N:DD6=DDDL2N:DY1=DYDL1N:DY6=DYDL2N:DPM1=DP 
MDL1N: DPM6=DPMDL2N : DYM1=DYMDL1N: DYM6=DYMDL2N : DRM1=D 
RMDL1N : DRM6=DRMDL2N 

IF DI(K) 5 0 THEN GOTO 1000 ELSE GOTO 1010 
DL1=DLDD1P: DL6=DLDD2P: DY2=DYDD1P: DY7=DYDD2P: DPM2=DP 
MDD1P: DPM7=DPMDD2P: DYM2=DYMDD1P: DYM7=DYMDD2P: DRM2=D 
RMDD1P: DRM7=DRMDD2P 

GOTO 1020 

DL1=DLDD1N: DL6=DLDD2N: DY2=DYDD1N: DY7=DYDD2N: DPM2=DP 
MDDIN : DPM7=DPMDD2N : DYM2=DYMDD1N : DYM7=DYMDD2N : DRM2=D 
RMDD1N : DRM7=DRMDD2N 

IF YL(K) > 0 THEN GOTO 1030 ELSE GOTO 1040 
DL2=DLDY1P: DL7=DLDY2P: DD2=DDDY1P: DD7=DDDY2P: DPM3=DP 
MDY1P: DPM8=DPMDY2P: DYM3=DYMDY1P: DYM8=DYMDY2P: DRM3=D 
RMDY1P: DRM8=DRMDY2P 

GOTO 1050 

DL2=DLDY1N: DL7=DLDY2N: DD2=DDDY1N: DD7=DDDY2N: DPM3=DP 
MDY1N: DPM8=DPMDY2N: DYM3=DYMDY 1N : DYM8=DYMDY2N : DRM3=D . 
RMDY 1N : DRM8=DRMDY2N 

IF PMI(K) > O THEN GOTO 1060 ELSE GOTO 1080 
DL3=DLDPM1P: DL8=DLDPM1P: DD3=DDDPM1P: DD8=DDDPM2P: DY3 
=DYDPM1P: DY8=DYDPM2P: DYM4=DYMDPM1P: DYM9=DYMDPM2P: DR 
M4=DRMDPM1P: DRM9=DRMDPM2P 

GOTO 1090 

DL3=DLDPM1N : DL8=DLDPM1N : DD3=DDDPM1N : DD8=DDDPM2N: DY3 
=DYDPM1N : DY8=DYDPM2N : DYM4=DYMDPM1N: DYM9=DYMDPM2N: DR 
M4=DRMDPM1N : DRM9=DRMDPM2N 

IF YM1(K) > 0 THEN GOTO 1100 ELSE GOTO 1110 
DL4=DLDYM1P: DL9=DLDYM2P: DD4=DDDYM1P: DD9=DDDYM2P: DY4 
=DYDYM1P: DY9=DYDYM2P: DPM4=DPMDYM1P: DPM9=DPMDYM2P: DR 
M5=DRMDYM1P: DRM10=DRMDYM2P 

Goro 1120 

DL4=DLDYM1N: DL9=DLDYM2N: DD4=DDDYM1N : DD9=DDDYM2N: DY4 
=DYDYM1N : DY9=DYDYM2N : DPM4=DPMDYM1N: DPM9=DPMDYM2N: DR 
M5=DRMDYM1N: DRM10=DRMDYM2N 

IF RM1(K) > 0 THEN GOTO 1130 ELSE GOTO 1140 
DL5=DLDRM1P: DL10=DLDRM2 P: DD5=DDDRM1P: DD10=DDDRM2P:D 
Y5=DYDRM1P:DY10=DYDRM2P: DPM5=DPMDRM1P: DPM10=DPMDRM2 
P: DYM5=DYMDRM1P: DYM10=DYMDRM2P 

GOTO 1150 
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1140 


1150 
1160 
Lo 
1220 
1230 
1250 
1260 
1270 
1320 
1550 
1340 
2330 
2339 


2340 
2350 
2360 


2370 


DL5=DLDRM1N:DL10=DLDRM2N:DD5=DDDRM1N:DD10=DDDRM2N:D 
Y5=DYDRM1N: DY10=DYDRM2N : DPM5=DPMDRM1N : DPM1LO=DPMDRM2 
N : DYM5=DYMDRMIN : DYM10=DYMDRM2N 


L2 (K)=L(K)-(DL1*D1 (K) )- (DL6* (D1 (K)*2) ) 
D2 (K) =D(K) -(DD1*L1(K) ) -(DD6*(D1(K) *2) ) 
'-==-COMPARE INTERACTION CALCULATIONS------ 
DFL=ABS (L1 (K) -L2 (K) ) : DFD=ABS (D1 (K) -D2 (K) ) 
IF DFL < .0005 AND DFD < .0005 THEN GOTO 1330 
IF DFL > .0005 THEN L1(K)=L2 (K) 
IF DFD > .0005 THEN D1(K)=D2 (K) 
GOTO 960 
L(K)=L2 (K) : D(K)=D2 (K) 
NEXT K 
'----PERFORM WEIGHT TARE CORRECTIONS-------- 
'--DRAG AND LIFT NEED TO BE CORRECTED FOR AXIS ORIENT 
ATION-- 
FOR K = 1 TO NOD% 
ALPHA = AOA(K)*(3.141593/180) 'CONVERT TO RADIANS 
DF (K) =(((-1*D(K) ) -(W*SIN (ALPHA) ))*COS(ALPHA)) + ((( 
=1*L(K) ) = (W*COS (ALPHA) ) ) *SIN (ALPHA) ) 
LF(K)=(((=2*L(K) ) = (W*COS (ALPHA) ) ) *COS (ALPHA) ) - ((( 
-1*L(K) )-(W*SIN (ALPHA) ) ) *SIN (ALPHA) ) 
NEXT K 
'=----SAVE REDUCED DATA TO FILE----------- 
CLS:PRINT "DATA REDUCTION COMPLETE, SAVING DATA TO "; 
FILES 
OPEN FILES FOR OUTPUT AS #2 
WRITE #2,NOD%,CONFIGS, CONF 
WRITE #2,0,W 
WRITE 42,ZL,ZD,ZV,ZPM,ZXVM,ZRM,ZAOA 
WRITE 42,CLL,CALD, CALY, CALPM, CALYM, CALRM, CALAOA 
FOR J = 1 TO NOD% 
WRITE 42,LF(J),DF(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) , AOA (J) 
NEXT J 


COMMON PDS, FDS,DX% 
CHAIN PDS+"MAIN.BAS",2110 
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260 


270 
280 
290 
300 
305 
310 
320 


330 
340 


350 
360 
IO 
380 
390 
400 
410 
415 
420 
430 


GOSUB 180 
COLOR 15,1,7: KEY OFF: CLS 
LOCATE 5,34: PRINT "*OPTIONS MENU*" 
PRINT 
PRINT 
PRINT TAB(30) "1. CREATE A NEW FILE" 
PRINT 
PRINT TAB(30) "2. EDIT EXISTING FILE" 
PRINT 
PRINT TAB(30) "3. INCREASE NUMBER OF TEST POINTS" 
PRINT TAB(30) " IN AN EXISTING FILE" 
PRINT 
PRINT TAB(30) "4. EXIT DATA EDITOR" 
LOCATE 20,5: INPUT "ENTER YOUR CHOICE (1,2,3,4)";REP% 
ON REP% GOSUB 390,280, 2530,480 
GOTO 30 


OETrONABASE 1: KEY OFF 

DIM LF(100),DF(100) ,YF(100),PM(100) ,RM(100) 
DIM YM(100) ,AOA(100) , COMMAND(10) 

BLANK2$ ü 

BLANK1S 
BLANKS Ë l 
COMMAND$ (1) — ''ef: COMMANDS(2) 
= "q" 

COMMAND$ (4) = "e": COMMAND$ (5) 
MUMJSGOMMANDS(7) = "q" 

RETURN 


"r": COMMAND$ (3) 


"s"; COMMANDS (6) = 


GLS? CÖLƏR 15,1,7 | 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME" ;FILES 
FILES = FDS+FILES 

GOSUB 2680 ' READ IN FILE 

PAGES = 1: COL% = l: SAVED% = 0: OLDXPOS% = 0: 
OLDYPOS% = 0 

GOSUB 1050 'PRINT BACKGROUND LINES 

ROWS = 3: GOSUB 1740 : GOSUB 1900 ‘HIGHLIGHT 
FIRST ROW AND COLUMN 

LOCATE 23,30: PRINT "LISTING": BEEP 

GOSUB 1430 'PRINT COLUMNS 1-7 AND Q VALUE 
LOCATE 23,20: PRINT BLANK2$ 

GOTO 470 


COLOR 15, 1,73 CIS 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME";FILE$ 

FILES = FD$ + FILES 

LOCATE 7,10: INPUT "ENTER TEST CONFIGURATION" ;CONFIG$ 
LOCATE 9,10: INPUT "ENTER THE NUMBER OF TEST POINTS 
PER RUN";NOD% 
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440 


450 
460 


470 
480 
481 
482 
490 
500 
S10 
520 
530 


540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 


680 
690 
700 
710 
720 
130 
740 
750 
760 
170 
780 


790 
800 


PAGES = Y: COLS = 1: SAVED% = 0: OLDXPOS% = 0: 
OLDYPOS$ = 0 

GOSUB 1050 "PRINT BACKGROUND LINES 

ROW% = 3: GOSUB 1740: GOSUB 1900 "HIGHLIGHT 
FIRST ROW AND COLUMN 

GOSUB 490 ' INSERT NEW DATA 

COMMON FDS,PD$,DX% 

CHAIN PDS+"MAIN.BAS",2110 


LOCATE 23,12: PRINT BLANKS 
LOCATE 23,12: INPUT REPS: IF REPS = "" THEN GOTO 530 
GOTO 570 
LOCATE 23,12: PRINT BLANKS: LOCATE 23,12: PRINT 
"INVALID COMMAND": BEEP 
FOR I = 1 TO 500 STEP 1 
NEXT I 
GOTO 500 
REP2$ = RIGHTS (REP$,1) 
Tl% = ASC(REP2$) 
IF T1% >= 43 AND T1% <= 57 THEN GOTO 650 
T1% = T1% OR 32: REP2$ = CHR$(T1%) 
GOSUB 2000 | 
IF VALID% = 0 THEN GOTO 500 ELSE VALID% = 0 
ON FLAG% GOSUB 1780,1670,2080,2340,2770,2170,2260 
IF REP2$ <> "e" THEN GOTO 500 ELSE GOTO 1040 
WHILE REP2$ <> "e" 
NEWDATA = VAL(REPS ) 
IF QUIT% = 1 OR OLDXPOS% = 0 THEN OUIT% = 0: GOTO 
700 ELSE 
LOCATE OLDXPOS%, OLDYPOS%, 0 
PRINT USING "####.##";OLDDATA 
LOCATE 23,12: PRINT BLANKS 
R% = ROW% - 2 
IF COL% = 1 THEN LF(R%+18*(PAGE%-1))= NEWDATA: 
YPOS% = 7: GOTO 790: ELSE 


IF COL% = 2 THEN DF(R%+18*(PAGE%-1) )=NEWDATA: YPOS% 


= 18: GOTO 790: ELSE 


IF COL% = 3 THEN YF(R%+18*(PAGE%-1) )=NEWDATA: YPOS% 
= 29: GOTO 790: ELSE 

IF COL% = 4 THEN PM(R%+18* (PAGE%-1) )=NEWDATA: YPOS% 
= 40: GOTO 790: ELSE 

IF COL% = 5 THEN RM(R%+18* (PAGE%-1) )=NEWDATA: YPOS% 
= 51: GOTO 790: ELS 

IF COL% = 6 THEN YM(R%+18*(PAGE%-1) )=NEWDATA: YPOS% 
= 62: GOTO 790: ELS 

IF COL% =7 THEN AOA(R%+18*(PAGE%-1))=NEWDATA: YPOS% 


= 71: GOTO 790: ELS 
XPOS% = ROW% 
LOCATE XPOS%,YPOS%: COLOR 0,10 
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810 
820 
830 


840 
850 
860 
870 
880 
890 
900 
910 
920 
930 


940 
950 
960 
970 
980 
990 
1000 
1010 


1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
BELO 
1120 
1550 
1140 
SO 
1160 
1670 
1150 
19 0 
1200 
1250 
1220 
1230 
1240 
1250 
1260 
1270 


BOCATE 1,53: 
LOCATET 1,62: 
LOCATE 1,74: 
LOCATE 21,1 


PRINT USING "####.##"; NEWDATA: MODIFIED% = 1 
SAVED% = 0: COLOR 15,1,7 

OLDYPOS% = YPOS%: OLDXPOS% = XPOS%: OLDDATA = 
NEWDATA 

FWD% = 1 

I% = ROW%+1: GOSUB 1690 'HIGHLIGHT NEW ROW 

IF I% > NOD%+2 THEN FWD% = 1 ELSE FWD% = O 

IF I% > 20 THEN FWD% = 1 ELSE FWD% = O 

GOSUB 1800 'HIGHLIGHT NEW ROW 

COLOR 15,1,7 

LOCATE 23,12: PRINT BLANKŠ 

LOCATE 23,12:INPUT REPS: IF REPS = "" THEN GOTO 930 
GOTO 950 

LOCATE 23,12: PRINT BLANKS: LOCATE 23,12: PRINT 
"INVALID COMMAND": BEEP 

GOTO 900 

REP2$ = RIGHTS (REPS,1) 

T1% = ASC(REP2S) 


IF Tl% >= 43 AND T1% <= 57 THEN GOTO 650 


T1% = T1% OR 32: REP2$ = CHRS(T1%) 
GOSUB 2000 "CHECK COMMANDS 

IF VALID% = 0 THEN GOTO 500 ELSE VALID% = 0 

ON FLAG% GOSUB 1780,1670,2080,2340,2770,2170,2260 
'COMMAND SUBROUTINES 

IF REP2$ <> "e" THEN GOTO 900 ELSE GOTO 1030 

WEND | 


RETURN 


COLOR 15,1,7 


LOCATE 2,1 

FORMES = 1 TO 80 : PRINT CHRS(220);: NEXT 
TOCATE 1,9: PRINT "LIFT" 

LOCATE 1,20: PRINT "DRAG" 

POCATE 1,32: PRINT "YAW" 

ESSATE” 1,41 "PBRINT' "PYTCA" 


PRINT "ROLL" 
PRINT "SAW M." 
PRINT "AOA" 


NEXT 


FOR I% = 1 TO 80 : PRINT CHR$(220);: NEXT 
FOR I% = 1 TO 22 

LOCATE I%,1: PRINT CHR$(222) 

LOCATE I%,80: PRINT CHRS(222) 
FOR I$ — 1 TO 21 


LOCATE I%,14: 
LOCATE I%,25: 
LOCATE 13,36: 
LOCATE I%,47: 


PRINT CHRŠ$ (179) 
PRINT CHRS$ (179) 
PRINT CHRŠ$ (179) 
PRINT CHR$ (179) 
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1280 
1290 
1300 
JO 
1:20 
1330 
1340 
1350 
1360 
1:70 


1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 


1460 
1470 
1480 


1490 
1500 
1510 


1520 
1:90 
1540 


1550 
1560 
1570 


1580 
15390 
1600 


LENO 
1620 
1630 


1640 
1650 
1660 
1670 
1680 
1690 


LOCATE 1$,58: PRINT CHRS(179) 
LOCATE I%,69: PRTNT CHRS(179) 


NEXT 

FOR I% = 1 TO 18 
LOCATE I%+2,2 
PRINT USING "444":1$+18*(PAGE%-1);: PRINT ":" 
NEXT 

COLOR 15,4 

LOCATE 22,15: PRINT "C COLUMN R ROW D DOWN E EXIT" 

PRINT " S SAVE U UP Q TUNNEL 


LOCATE 22,46: 
SPEED" — 
COLOR 14,0: 
COLOR 15,1,7 
LOCATE 22,2: 


LOCATE 23,44: PRINT CONFIGS 


PRINT "10:27": UZCATE 22712. 2EEELE””7” 


LOCATE 23,3: PRINT "COMMAND :" 

RETURN 

.--------- PRINT COLUMNS 1-7 AND Q VALUE--------- 

FOR I% = 1 TO 18 
LOCATE I%+2,7: PRINT USING "###44.4#4";LF(I%+18* (PAG 
E%-1)) 
NEXT 

FOR I% = 1 TO 18 | 
LOCATE I%+2,18: PRINT USING "####.44";DF(I%+18*(PA 
GES-1)) 
NEXT 


FOR I% = 1 TO 18 
LOCATE 1%+2,29: 
GES-1)) 

NEXT 

FOR 1% = 1 TO 18 
LOCATE 1%+2,40: 
GE%-1) ) 

NEXT 

FOR 13 = 1 TO 18 
LOCATE I%+2,51: 
GES-1)) 

NEXT 

FOR I% = 1 TO 18 
LOCATE I%+2,62: 
GE%-1) ) 

NEXT 

FOR I% = 1 TO 18 
LOCATE 1%+2,71: 
AGE%-1)) 

NEXT 
LOCATE 22,7:PRINT USING "###";Q 


PRINT USING "####.##" 7; YF (I%+18* (PA 


PRINT USING "4444. 44" ;PM(I$+18* (PA 
PRINT USING "4ätd. 44" ;RM(1$+18* (PA 
mid HUW VM(IS+18* (PA 


PRINT USING 


PRINT USING "####.##";AOA(I%+18*(P 


FWD% = VAL(REP$): IF FWD% = 0 THEN FWD% = 1 
OLDROWS = ROWS: ROW% = ROW% + FWD% 
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1700 
1710 
1720 
1730 


1740 
1759 
1760 
1770 
1780 
17290 
1800 
1810 
1820 
18.30 


1840 


1850 


1860 


1870 


1880 


1590 


1900 
10 


1920 
1930 
1940 
1850 
1960 
1970 
1980 
1990 
2000 
2010 
2020 


2030 
2040 


IF ROW% > NOD%+2 THEN ROW% = 3 

IF ROW% > 20 THEN ROW% = 3 

IF ROW% < 3 THEN ROW% = 3 

LOCATE OLDROW%,2: PRINT USING "###" ;OLDROW%-2+18* (PAG 


ES-1);: PRINT Ws 
COLOR 14,4: LOCATE ROWZ,2 


PRINT USING "444";ROW%-2+18*(PAGE$-1);: 
COLOR 15,1,7 


PRINT ":" 


RETURN 
Ur COLUMN INDEXING----------------- 
FWD% — VAL(REP$): IF FWD% = O THEN FWD% = 1 
OLDCOL% = COL%: COL% = COL% + FWD% 
IF COL% > 7 THEN COL% = 1 
IF COL% < 1 THEN COL% = 7 
IF OLDCOL$ = 1 THEN LOCATE 1,9: PRINT "LIFT": GOTO 
1900 ELSE 
IF OLDCOL% = 2 THEN LOCATE 1,20: PRINT "DRAG": GOTO 
1900 ELSE 
IF OLDCOL% = 3 THEN LOCATE 1,32: PRINT "YAW": GOTO 
1900 ELSE 
IF OLDCOL% = 4 THEN LOCATE 1,41: PRINT "PITCH": GOTO 
1900 ELSE 
IF OLDCOL% = 5 THEN LOCATE 1,53: PRINT "ROLL": GOTO 
1900 ELSE 
IF OLDCOL% = 6 THEN LOCATE 1,62: PRINT "YAW M.": GOTO 
1900 ELSE 
IF OLDCOL$ = 7 THEN LOCATE 1,74: PRINT "AOA": GOTO 
1900 ELSE 
COLOR 14,4 
IF COL% = 1 THEN LOCATE 1,9: PRINT "LIFT": GOTO 1980 
ELSE 
IF COL% = 2 THEN LOCATE 1,20: PRINT "DRAG": GOTO 1980 
ELSE | 
IF COL% = 3 THEN LOCATE 1,32:: PRINT "YAW": GOTO 1980 
ELSE | 
IF COL% = 4 THEN LOCATE 1,41: PRINT "PITCH": GOTO 1980 
ELSE 
IF COL% = 5 THEN LOCATE 1,53: PRINT "ROLL": GOTO 1980 
ELSE 
IF COL% = 6 THEN LOCATE 1,62: PRINT "YAW M.": GOTO 
1980 ELSE 
IF COL% = 7 THEN LOCATE 1,74: PRINT "AOA": GOTO 1980 
ELSE 
COLOR 15,1,7 
RETURN 
.------ CHECK COMMANDS--------------------- 
I% = 0 
WHILE I% < 7 AND VALID% = 0 

I% = I% + 1 

IF REP2$ = COMMAND$(I%) THEN FLAG% = 1%: VALID% = 1 
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2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 


2140 
2150 
2160 
2170 
2180 
2190 
2200 
2257 
2220 


2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2610 
2520 


WEND 

IF VALID$ = 0 THEN GOSUB 2420 "ERROR MESSAGE 
RETURN 
! ーーーーーーーーーー SCROLL DOVN----------------- 
FWD$ = VAL(REP$): IF FWD% = O THEN FWD% = 1 
PAGE% = PAGE% + FWD% 


IF PAGES < 1 THEN PAGES = 1 


GOSUB 1050 "PRINT BACKGROUND LINES 

ROWS = 3: COL% = 1: GOSUB 1740: GOSUB 1900 'HIGHLIGHT 
FIRST ROW AND COLUMN 

GOSUB 1430 'PRINT COLUMNS 1-7 AND Q VALUE 

OLDXPOS% = 0 | 

RETURN 

.--------- SCROLL UP--------------------- 

FWD% = VAL(REPS): IF FWD% = 0 THEN FWD% = 1 

PAGES = PAGE% - FWD% 

IF PAGES < 1 THEN PAGES = 1 

GOSUB 1050 'PRINT BACKGROUND LINES 

ROW% = 3: COL% = 1: GOSUB 1740: GOSUB 1900 'HIGHLIGHT 
FIRST ROW AND COLUMN 

GOSUB 1430 ' PRINT COLUMNS 1-7 AND Q VALUE 

OLDXPOS% = 0 


-------- PRINT OUT Q VALUE--------------- 
LOCATE 23,12: PRINT BLANKS 

LOCATE 23,12: INPUT "Q =":Q 

LOCATE 23,12: PRINT BLANKS 

corso 10 
LOCATE 22 7; 
COLOR 15,1 7 


PRINT USING "###";Q 


IF SAVED% = 1 OR MODIFIED% = 0 THEN GOTO 2410 
LOCATE 23,20: PRINT BLANK2$: BEEP ` 
LOCATE 23,20: INPUT "SAVE FILE (Y/N)";REPS 


IF REPS = "Ы" OR REPS = "n" THEN GOTO 2410 
GOSUB 2770 "SAVE FILE 
LOCATE 23,20: PRINT BLANK2$ 


LOCATE 23,20: PRINT BLANK2$ 


LOCATE 23,20: PRINT " INVALID COMMAND" 
FOR F = 300 TO 500 STEP 100 
SOUND F,2 
SOUND 32767,2 
NEXT 
FOR I = 1 TO 500 STEP 1 
NEXT I 
LOCATE 23,20: PRINT BLANK2$ 
RETURN 


Figure A.12 ADATA.BAS (cont.) 


64 


2540 
2550 


2560 
2570 
2580 
2590 
2610 


2020 


2630 
2660 


2670 
2680 
2690 
2700 
2710 
2714 
2715 
2720 
2730. 
2740 
2750 
2760 
2770 
2780 


2790 
2800 
2810 
2814 
2915 
2820 
2830 
2840 
2850 
2860 


. .— —-— ADD TESIS POTNTS TO PILEE=--------- 
CLS 
LOCATE 3,5: PRINT"THIS OPTION ALLOWS YOU TO INCREASE 
THE NUMBER OF TEST" 
LOCATE 4,5: PRINT"POINTS OF AN EXISTING FILE" 
LOCATE 7,5: INPUT"VHAT IS THE NAME OF YOUR FILE"”FILEŞ 
FILEŞ — FDŞ + FILES 
GOSUB 2680 
LOCATE 10,5: INPUT"WHAT IS THE NEW NUMBER OF TUNNEL 
RUNS" ;NOQ% 
LOCATE 13,5: INPUT"WHAT IS THE NEW NUMBER OF TEST 
POINTS PER RUN" ;NOD% 
GOSUB 2770 
LOCATE 20,5: PRINT"NOW YOU CAN CALL UP YOUR FILE AND 
ADD IN THE NEW POINTS" 
RETURN 


OPEN FILES FOR INPUT AS +1 

INPUT #1,NOD%, CONFIGS, CONF 

INPUT #1,9,W 

INPUT £1,ZL,ZD,ZY,ZPM,ZYM,ZRM,ZAOA 

INPUT #1,CLL,CALD, CALY, CALPM, CALYM, CALRM, CALAOA 
FOR J = 1 TO NOD% 
INPUT 41,LF(J) ,DF(J) ,YF(J) ,PM(J) ,RM(J) ,YM(J) ,AOA(J) 
NEXT J 

CLOSE #1 

RETURN 

:.------------ SAVE FILE-------------------- 

LOCATE 23,20: PRINT BLANK2$: LOCATE 23/30: PRINT 

"SAVING FILE" 

OPEN FILES FOR OUTPUT AS #2 

WRITE #2,NOD%, CONFIGS, CONF 

WRITE £2,0,W 

WRITE +2,ZL,ZD,ZY, ZPM, ZYM, ZRM, ZAOA 

WRITE #2,CLL,CALD, CALY, CALPM, CALYM, CALRM, CALAOA 

FOR J = 1 TO NOD% 
WRITE #2,LF(J),DF(J),YF(J),PM(J) ,RM(J) ,YM(J) ,AOA(J) 
NEXT J l 

CLOSE 42 

RETURN 
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REM----COMP.BAS; CALCULATE Cd, Cl AND E. F. A.----- 

REM 

DIM LF(100),DF(100),YF(100),PM(100) ,YM(100),RM(100), 

AOA(100) 

DIM CD(100) ,CL(100) , EFA(100) 

COLOR 15,1,4:KEY OFF: CLS 

PRINT "WHAT IS THE NAME OF YOUR FILE CONTAINING THE"; 
PRINT " CALIBRATED" 

INPUT "WIND TUNNEL DATA";FILES 

FILES = FD$+ FILES 


GOSUB 680 'READ IN. DATA FROM FILE 
REM -- FOR THE FOLLOWING, AREAL IS MEASURED IN SQUARE 
FEET 


AREAL =60/144 
FOR J = I TO NOD% 
Chita LF (J) / (Q*AREAL) 
CD (J) DF(J)/(Q*AREAL) 
EFA(J) =LF(J)/Q 
NEXT J 
(== = s. STORE DATA FOR Cam S. Cl a 
INPUT" WHAT IS THE NAME FOR THE FILE TO STORE Cad VS. 
Cl DATA";NS 
Ys="cq" * Хӱ="С1 " 
OPEN FD$+N$ FOR OUTPUT AS #1. 
WRITE #1,NOD%,Y$,X$, CONFIGS, CONF 


WRITE 41,0 
FOR J = 1 TO NOD% 
WRITE #1,CD(J) ,CL(J) 
NEXT J 


----------- STORE DATA FOR Cd VS. Cl“C1------------ 
INPUT"WHAT IS THE NAME FOR THE FILE TO STORE Cd VS. 
Ci?) DATA ANS 
yss" са" 5 = 0 
OPEN FD$+N$ FOR OUTPUT AS #2 
WRITE #2,NOD%,Y$,X$,CONFIGS , CONF 
WRITE 42,0 

FOR J = 1 TO NOD% 
B= CL(I)*CL(J) 
WRITE #2,CD(J),B 
NEXT J 


INPUT"WHAT IS THE NAME FOR THE FILE TO STORE THE Cl 
VS. AOA DATA";N$ 
YSZEGIN = VAIDA 
OPEN FDS+N$ FOR OUTPUT AS +2 
WRITE %2,NOD%,YS,XS$,CONFIGS,CONF 
VRITE 42,0 
FOR J = 1 TO NODŠ 
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530 WRITE #2,CL(J) ,AOA(J) 


540 NEXT J 
550 CLOSE #2 
560 ge = STORE DATA FOR E.F.A. VS AOA----------- 


561  INPUT'WHAT IS THE NAME FOR THE FILE TO STORE THE 
.F.A. VS. AOA DATA"?N$ 

562 Y$-"E.F.A.":X$-"AOA" 

570 OPEN FDS+N$ FOR OUTPUT AS #1 

580 WRITE #1,NOD%,Y$,X$, CONFIGS, CONF 

590 WRITE 41,0 


600 FOR J = 1 TO NOD% 

610 WRITE #1,EFA(J) ,AOA(J) 

620 NEXT J 

630 CLOSE 41 

640 (-------- RETURN TO MAIN PROGRAM--------------------- 


650 COMMON FDS, PDS, Dx% 

660 CHAIN PDS+"MAIN.BAS",2110 

670 (------------ READ IN DATA FROM CALIBRATED FILE----- 
680 OPEN FILES FOR INPUT AS #2 

690 INPUT %2,NOD%, CONFIGS$,CONF 

700 INPUT 42,0,W 

701 INPUT 42,ZL,ZD,ZV,ZPM,ZVM,ZRM,ZAOA 

702 INPUT 42,CLL,CALD,CALV,CALPM, CALVM, CALRM, CALAOA 


710 FOR J = 1 TO NOD% 

720 INPUT #2,LF(J),DF(J),YF(J),PM(J),RM(J),YM(J), 
AOA ( け ) 

730 NEXT J 


740 CLOSE #2 
750 RETURN 
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GOSUB 180 
COLOR 15,1,7: KEY OFF: CLS 
LOCATE 5,34: PRINT "*OPTIONS MENU*" 
PRINT 
PRINT 
PRINT TAB(30) "1. CREATE A NEW FILE" 
PRINT 
PRINT TAB(30) "2. EDIT EXISTING FILE" 
PRINT 
PRINT TAB(30) "3. INCREASE NUMBER OF TEST POINTS" 
PRINT TAB(30) " IN AN EXISTING FILE" 
PRINT 
PRINT TAB(30) "4. EXIT DATA EDITOR" 
LOCATE 20,5: INPUT "ENTER YOUR CHOICE (1,2,3,4)";REP% 
ON REP% GOSUB 380,270,2200,490 
GOTO 30 


OPTION BASE 1: KEY OFF 

DIM Y(100) ,X(100) , COMMAND(10) 
BLANK2$ 
BLANK1$ 
BLANKS 
COMMANDS (1) = "c":COMMANDS (2) = "r": COMMAND$(3) = "a" 
COMMAND$ (4) = "e": COMMANDS(5) = "s": COMMANDS (6) = 
"üns COMMANDŞ (7) — üçü. 

RETURN 


CES: COFORSI5/1,/7 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME"; FILES 
FEBES = FDS +" FILES 

GOSUB 2350 ' READ IN FILE 

PAGES = 1: COL% = 1l: SAVED% = 0: OLDXPOS% = O: 
OLDYPOS% = 0 

GOSUB 1010 'PRINT BACKGROUND LINES 

ROW% = 3: GOSUB 1520 : GOSUB 1620 ‘HIGHLIGHT 
FIRST ROW AND COLUMN 

LOCATE 23,30: PRINT ETSTENG BEEP 

GOSUB 1360 'PRINT COLUMNS l-7 AND Q VALUE 
LOCATE 23,20: PRINT BLANK2$ 

GOTO 480 


COLOR 15,1,7: CLS 

LOCATE 5,10: INPUT "ENTER YOUR FILE NAME";FILES 

FILRS = FDS + FILES 

LOCATE 7,10: INPUT "ENTER TEST CONFIGURATION" ; CONFIGS 
LOCATE 9,10: INPUT "ENTER THE NUMBER OF TEST POINTS 
PER RUN" ;NOD% 

LOCATE 11,10: INPUT "ENTER TITLE FOR X-VALUES (i.e. 
CL, AOA)"”X1$ 
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LOCATE 13,10: INPUT "ENTER TITLE FOR Y-VALUES (i.e. 
ca merj": Y1S 
PAGES = 1: COL% = 1: SAVED% = 0: OLDXPOS% = O: 
OLDYPOS% = 0 
GOSUB 1010 ‘PRINT BACKGROUND LINES 
ROW% = 3: GOSUB 1520: GOSUB 1620 ‘HIGHLIGHT FIRST 
ROW AND COLUMN 
GOSUB 500 ' INSERT NEW DATA 
COMMON FDS, PDS$,DxX% 
CHAIN PD$+"MAIN.BAS",2110 
RETURN 
——— oo CHANGE OR INSERT DATA--------------- 
LOCATE 23,12: PRINT BLANK$ 
LOCATE 23,12: INPUT REPS: IF REPS = "" THEN GOTO 540 
GOTO 580 
LOCATE 23,12: PRINT BLANKS: LOCATE 23,12: PRINT 
"INVALID COMMAND": BEEP 
FOR I = 1 TO 500 STEP 1 
NEXT I 
GOTO 510 
REP2$ = RIGHTS(REPS,1) 
T1% = ASC(REP2$) 
IF T1% >= 43 AND T1% <= 57 THEN GOTO 660 
T1% = T1% OR 32: REP2S = CHRS(T1%) 
GOSUB 1670 | 
IF VALID% = 0 THEN GOTO 510 ELSE VALID% = 0 
ON FLAG% GOSUB 1550,1450,1750,2010,2440,1840,1930 
IF REP2S <> "e" THEN GOTO 510 ELSE GOTO 1000 
WHILE REP2S <> "et 
NEWDATA = VAL(REPS) 
EUI. =I OR OLDXPOS% = O THEN QUITS = "0: GOTO 
710 ELSE 
LOCATE OLDXPOS%,OLDYPOS%,0 
PRINT USING "####.##";OLDDATA 
LOCATE 23,12: PRINT BLANKS 
R% = ROW% - 2 | 
IF COL% = 1 THEN Y(R%+18*(PAGE%-1) )=NEWDATA: YPOS% 
= 18: GOTO 750: ELS 
IF COL = 2 THEN X(R%+18*(PAGE%-1) )=NEWDATA: YPOS$% 
= 43: GOTO 750: ELSE 
XPOS% = ROW% 
LOCATE XPOS%,YPOS%: COLOR 0,10 
PRINT USING "####.##"; NEWDATA: MODIFIED% = 1 
SAVED% = 0: COLOR 15,1 
OLDYPOS% = YPOS%: OLDXPOS% = XPOS%: OLDDATA = 
NEWDATA 
FWD% = 1 
I% = ROW%+1: GOSUB 1470 "HIGHLIGHT NEW ROW 
IF I% > NOD%+2 THEN FWD% = 1 ELSE FWD% = 0 
IF I% > 20 THEN FWD% = 1 ELSE РИОЗ = О 
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GOSUB 1570 
COLOR 15,1,7 
LOCATE 23,12: 


"HIGHLIGHT NEV ROV 


PRINT BLANKŞ 


LOCATE 23,12: INPUT REP$: IF REPS = "" THEN GOTO 
890 

GOTO 910 

LOCATE 23,12: PRINT BLANKS: LOCATE 23,12: PRINT 
"INVALID COMMAND": BEEP 

GOTO 860 

REP2$ = RIGHTS (REPS, 1) 


T1% = ASC(REP2$) 


IF Tl% >= 43 AND T1% <= 57 THEN GOTO 660 


T1% = T1% OR 32: REP2$ = СНЕ5(Т15) 
GOSUB 1670 'CHECK COMMANDS 
IF VALID% = 0 THEN GOTO 510 ELSE VALID% = 0 


ON FLAG% GOSUB 1550,1450,1750,2010,2440,1840,1930 
' COMMAND SUBROUTINES 
IF REP2$ <> "e" THEN GOTO 860 ELSE GOTO 990 


WEND 

RETURN 

— ニニ ニニ ーPRTNT BAGKGROUNEMERN S Sa a = 
COLOR 15,1,7 

CLS 
LOCATE 2,1 
FOR 1% = 1 TO 80 : PRINT CHR$(220);: NEXT 
LOCATE 1,18: PRINT "Y = ";Y1$ 

LOCATE 1,439 PRINT VX = 515 

LOCATE 1,68: PRINT "COMMANDS" 

LOCATE 21 
FOR I% = 1 TO 80 : PRINT CHR$(220);: NEXT 
FOR I% = 1 TO 22 


LOCATE I%,1: 
LOCATE 1%,80: 
NEXT 
FOR 1% = 1 TO 21 
LOCATE I%,28: 
LOCATE I%,53: 
NEXT 
FOR 1% = 1 TO 18 
LOCATE 1%+2,2 


PRINT CHRS (222) 
PRINT CHR$ (222) 


PRINT CHR$ (179) 
PRINT CHR$ (179) 


PRINT USING "ĘĘĘ";1%+18*(PAGE%-1);: PRINT ":" 
NEXT 
COLOR 15,4 
LOCATE 3,66: PRINT "C COLUMN " 
LOCATE 4,66: PRINT "R ROW i 
LOCATE 5,66: PRINT "D DOWN 4 
LOCATE 6,66: PRINT "E EXIT 
LOCATE 7,66: PRINT "S SAVE U 
LOCATE 8,66: PRINT ''U UP " 
LOCATE 9,66: PRINT "Q TUNNEL SPEED" 
Figure A.14 BDATA.BAS (cont.) 
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1310 COLOR 14,0: LOCATE 23,44: PRINT CONFIGS 
ern, COLOR LS,1,7 
1320 LBOCANE 282: PRUENT "10 ="; LOCATE 22,12: PRINT "|" 
1340 LOCATE 23,3: PRINT "COMMAND :" 
1360 RETURN 
' 


1360 "--------- PRINT COLUMNS 1-7 AND Q VALUE--------- 

1200 BƏR I£ = 1 TO 18 

1380 LOCATE I%+2,18: PRINT USING "####.##"7Y(I13+18* (PAG 
E%-1)) 

1390 NEXT 

00 FOR I% = 1 TO 18 

1410 LOCATE I%$+2,43: PRINT USING "####.##";X(1I%+18* (PAG 
E%-1)) 

1420 NEXT 


1430 LOCATE 22,7:PRINT USING "###"7@ 

1440 RETURN 

1450 "---------- ROV INDEXING---------------- 

1460° FWD% = VAL(REPS): IF FWD% THEN FWD% = 1 

1470 OLDROW% = ROW%: ROW% = ROW F 

1480 IF ROW% > NOD%+2 THEN ROW% = 3 

1490 IF ROW% > 20 THEN ROW% = 3 

1500 IF ROW% < 3 THEN ROW% = 3 | 

1510 LOCATE OLDROW%,2: PRINT USING "###" ;OLDROW%-2+18* (PAG 
Eżejdw : RINT "ес" 

1520 COLOR 14,4: LOCATE ROW%,2: PRINT USING "###" ;ROW%-2+1 
8*(PAGE%-1) ;:PRINT ":" 

1530 COLOR 15,1,7 

1540 RETURN 

1550 "---------- COLUMN INDEXING----------------- 


oo İl 
+ O 


1560 FWD% = VAL(REPS): IF FWD% = O THEN FWD% = 1 
1570 OLDCOL% = COL%: COL% = COL% + FWD% 


1580 IF COL% > 2 THEN COL% = 1 
1590 IF COL% < 1 THEN COL% = 2 


1600 IF OLDCOL% = 1 THEN LOCATE 1,18: PRINT "Y = ";yY1S: 
GOTO 1620 ELSE 

1610 IF OLDCOL% = 2 THEN LOCATE 1,43: PRINT "X = ";X1$: 
GOTO 1620 ELSE 

1520 COLOR 14,4 | 

1630 IF COL$ = 1 THEN LOCATE 1,18: PRINT "Y = ";Y1$: GOTO 
1650 ELSE 

1640 IF COL$ = 2 THEN LOCATE 1,43: PRINT "X = ";X1$: GOTO 
1650 ELSE 

1650 COLOR 15,1,7 

1660 RETURN 

1670 "------- CHECK COMMANDS--------------------- 

1680 I% = 0 

1690 WHILE 19 “= 7 AND VALID% = 0 

1700 I% = I% + 1 

1710 IF REP2$ = COMMANDS(I%) THEN FLAGS = 1%: VALID% = 1 

1720 WEND 
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IF VALID% = 0 THEN GOSUB 2090 "ERROR MESSAGE 


RETURN 

! ーーーーーーーーーー SCROLL DOVN----------------- 
FWD$ = VAL(REPS): IF FWD$ = 0 THEN FWD% = 1 
PAGE% = PAGE% + FWD% 


IF PAGE? € 1 THEN PAGES = 1 


GOSUB 1010 'PRINT BACKGROUND LINES 

ROWS = 3: COL% = 1: GOSUB 1520: GOSUB 1620 'HIGHLIGHT 
FIRST ROW AND COLUMN 

GOSUB 1360 'PRINT COLUMNS 1-7 AND Q VALUE 

OLDXPOS% = O 

RETURN 

"ーー SCROLL UP——— ee 

FWD = VAL(REPS): IF FWD% = 0 THEN FWD% = 1 

PAGES = PAGES - FWD% 

IF PAGES < 1 THEN PAGES = 1 

GOSUB 1010 'PRINT BACKGROUND LINES 

ROW る = 3: COL% = 1: GOSUB 1520: GOSUB 1620 'HIGHLIGHT 
FIRST ROW AND COLUMN 

GOSUB 1360 ' PRINT COLUMNS 1-7 AND Q VALUE 

OLDXPOS% = 0 


ーーーーーーーー PRINT OUT Q VALUE--------------- 
LOCATE 23,12: PRINT BLANKS 

LOCATE 23,12: INPUT 1O =";0 

LOCATE 23,12: PRINT BLANK$ 

COLOR 0,10 f 

LOCATE 22,7: PRINT USING "###";Q 

COHOR IS 1 7 


O THEN GOTO 2080 
PRINT BLANE2$: BEEP 
INPUT "SAVE FILE (Y/N)";REPS 


IF SAVED% = 1 OR MODIFIED% = 
LOCATE 23,20: 
LOCATE 23,20: 


IF REPS = "N" OR REPS = "n" THEN GOTO 2080 
GOSUB 2440 'SAVE FILE 
LOCATE 23,20: PRINT BLANK2$ 


LOCATE 23,20: PRINT BLANK2$ 


LOCATE 23,20: PRINT " INVALID COMMAND" 
FOR F = 300 TO 500 STEP 100 
SOUND F,2 
SOUND 32767,2 
NEXT 
FOR I = 1 TO 500 STEP 1 
NEXT I 
LOCATE 23,20: PRINT BLANK2$ 
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CLS 

LOCATE 3,5: PRINT"THIS OPTION ALLOVS YOU TO INCREASE 
THE NUMBER OF TEST" 

LOCATE 4,5: PRINT"POINTS OF AN EXISTING FILE" 


LOCATE 7,5: INPUT"WHAT IS THE NAME OF YOUR FILE";FILE$ 
FILES = FDS + FILES 

GOSUB 2350 

LOCATE 13,5: INPUT"WHAT IS THE NEW NUMBER OF TEST 
POINTS PER RUN" ;NOD% 

GOSUB 2440 

LOCATE 20,5: PRINT"NOW YOU CAN CALL UP YOUR FILE AND 
ADD IN THE NEW POINTS" 


OPEN FILES FOR INPUT AS #1 
INPUT #1,NOD%,Y1$,X1$, CONFIGS, CONF 
INPUT 41,0 
FOR J = 1 TO NOD% 
INPUT 41,Y(J),X (9) 
NEXT J 
CLOSE +1 


LOCATE 23,20: 
"SAVING FILE" 
OPEN FILES FOR OUTPUT AS #2 
WRITE #2,NOD%,Y1$,X1$,CONFIGS , CONF 
WRITE 42,0 
FoR J = 1 TO NODS 

WRITE +2,Y(J),X(J) 

NEXT J 
CLOSE #2 
RETURN 


PRINT BLANK2$: LOCATE 23,30: PRINT 
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! ーーーーーーーーーーー 一 PLOTTING ROUTINE------------------ 
DIM X(300) ,Y(300) ,B(300) ,C(300) ,D(300) 
KEY OFF: SCREEN 2 
HDPOS% = 44: B$ = " " 
GOSUB 260 
OPTS = CURVES(1) 
ERASE X,Y,FILESTKS, CURVES,B,C,D 
CHAIN PDS+"MAIN. BAS" ,2110,ALL 
-——--------- READ IN PILE--- 1 Je E 
ən FILES FOR INPUT AS #2 
INPUT 42,NOD%,Y1$,X1$,H$,CONF 
INPUT +2,0 
FOR I% = 1 TO NOD% 
INPUT #2 7717). Xx (15) 
NEXT 
CLOSE +2 
RETURN 


FOR L% = 1 TO NOF% 
FILES = FILESTKS(L%) 
GOSUB 160 
IF L%=1 THEN XMIN=X(1): XMAX=XMIN: YMIN=Y (1): 
YMAX=YMIN 
FOR I%=1 TO NOD%3 
IF X(I%)<XMIN THEN XMIN=X(1%) ELSE IF X(1%)> 
XMAX THEN XMAX=X (18) 
IF Y(I%)<YMIN THEN YMIN=Y(1%) ELSE IF Y(1%)> 
YMAX THEN YMAX=Y (1%) 


NEXT 13 
NEXT L% 


CLS 

PRINT "AFTER CURVES ARE PLOTTED YOU HAVE THE OPTION 
OF OBTAINING A HARD" 

PRINT “COPY BY PRESSING CTKEZPBİ 7700 CURVE MAD 
REMAIN ON THE SCREEN" 

PRINT eit YOU PRESS THE ForkEX 


PRINT 

INPUT "HEADDING FOR PLOT =":;HEADS 

T1%=LEN(HEADS): IF T1%>60 THEN HEADS=LEFTS(HEADS, 60) 
INPUT "X TITLE FOR THE PLOT =":XTITLES 


TIlSSLEN(XTITLES): IF T15>60 THEU re HEEO—EBETO(XULDUR 
9,60) 

INPUT YY TITLE FOR THE PLOT =" YTI LELE; 
T1%=PEN(YTITLES): IF Tis>60 THEN YT EEE unu. 
9,60) 


PRINT "MINIMUM X VALUE = (DEFAULT:";XMIN;"]";: INPUT 
T1S 

IF Tl$ <> "" THEN TEMP=VAL(T1$) ELSE TEMP=XMIN 

IF TEMP<=XMIN THEN OK%=1 ELSE OK%= 

IF Tl$ <> "" THEN TEMP3 =XMIN: XMIN = TEMP 
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PRINT "MAXIMUM X VALUE = (DEFAULT:";XMAX;"]";: INPUT 
PES 

IF TIST <> "un XMAX-VAL(TIS) 

PRINT "MINIMUM Y VALUE = [DEFAULT:";YMIN;"]";: INPUT 
TIS 

IF T1$ <> "" THEN YMIN=VAL(T1$) 

PRINT "MAXIMUM Y VALUE = [DEFAULT:";YMAX;"]";: INPUT 
TIS 

IF T1$ <> "" THEN VMAX-VAL(TIS) 


RANGE=XMAX-XMIN 

IF OK%=1 THEN PTR%=1: STKPTR%=1: PGSTK%(1)=1 ELSE 
GOSUB 580 | 

CLS: SCREEN 2 


GOSUB 840 

RETURN 

! ーーーーーーーーーー SET PAGES----------------------------- 
TEMP2=(TEMP-TEMP3) /RANGE+1: PG%=INT(TEMP2) 


STKPTR%=PG%+1 
FOR 1% = 1 TO STKPTR% 
T2-TEMP-RANGEX (Is-1) 
J%=1 
WHILE X(J%) <T2 
J%=J%+1 
WEND 
PGSTK% (PG%-1%+2)=J% 
NEXT I$ 
PTR%=PGSTK% (STKPTR%) 
RETURN 


LINE (ORIGINX% , ORIGINY%) - (ORIGINX$+7* INTERVALX$% , TEMP4% 
) LB 
X%=ORIGINX% 
FOR I% = 1 TO 8 
LINE (X%,162)-(X%,165) 
X%=X$+INTERVALX% 
NEXT 
T1%=ORIGINX%-8: T2%=ORIGINY3S 
FOR I% = 1 TO 5 
LINE (T1%,T2%)-(ORIGINX%,T2%) 


T2%=T2%+INTERVALY% 

NEXT 
RETURN 
Ca ss = PECOT TAA ES = ——o——————-—————— 
ORIGINY%=12: ORIGINX%=72: INTERVALY%=30: INTERVALX%=80 


TEMP 4 %=ORIGINY$+5* INTERVALY% 
VIEW(0,0)-(639,163): CLS: VIEW 
GOSUB 710 

IF LEFTS(GRD$,1)="G" THEN GOSUB 2400 
X%=ORIGINX%: Y%=ORIGINY3 

FOR L%=1 TO NOF% 
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IF NOF%=1 THEN GOTO 960 
FILES-FILESTKS (L$) 
GOSUB 160 
XINDEX-XMIN 
GOSUB 1410 
IF L%=1 THEN GOSUB 1120 
IF L$ <> NOF$ THEN PTR%=PGSTK% (STKPTRS) 
NEXT L% 
ON KEY(1) GOSUB 1006 
ON KEY(2) GOSUB 1009 
KEY(2) OFF: KEY(1) ON: KEY(2) ON 
IF FLAG = 1 THEN GOTO 1100 
GOTO 1001 
LOCATE 23,1: PRINT "HI" 
RETURN 
FLAG = 1 
RETURN 
FLAG = 0 


--------- PLOT ZERO LINE------------------ 

IF (YMAX*YMIN <0) THEN TEMP=(-1*YMIN/ (YMAX-YMIN) ) *150 
ELSE GOTO 1150 

DELTAY%=CINT (TEMP) 

LINE (ORIGINX%,TEMP4%-DELTAY%) -(ORIGINX%+7*INTERVALX$%, 
TEMP4%-DELTAY3) 

IF (XMAX*XMIN <0) THEN TEMP=(-1*XMIN/ (XMAX-XMIN) ) *550 
ELSE GOTO 1190 

DELTAX%=CINT (TEMP) 

LINE (ORIGINX%+DELTAX% , ORIGINY% ) - (ORIGINX$+DELTAX$ , TEM 


YINDEX=YMAX 
TEMP = (YMAX-YMIN) /5 
P%=LEN (HEADS): P%=HDPOS%-(P%*.5) 
LOCATE 1,1: PRINT TAB(P%) HEADS 
FOR I% = 1 TO 6 
PRINT USING "##4#4.444";YINDEX 
IF (I$= 3) THEN FOR J% — 1 TO 2: PRINT: NEXT J$: 
GOTO 1260 
IF (18<>6) THEN FOR J% = 1 TO 3: PRINT: NEXT J% 
YINDEX=YINDEX-TEMP 
NEXT I% 
XINDEX = XMIN: XINCR=(XMAX-XMIN) /7 
LOCATE 22,1 
PRINT USING "######44.4">XINDEX; 
FOR 1% = 1 TO 7 
XINDEX = XINDEX + XINCR 
PRINT USING "######444.#";XINDEX; : NEXT I3 
P%=LEN(XTITLES): P%=HDPOS%-(P%*.5) 
LOCATE 23,1: PRINT TAB(P%) XTITLES 
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T1%=LEN(YTITLES): P%=11-(T1%*.5) 

FOR 1% = 1 TO Tl% 
T1$=MIDS(YTITLE$,I%,1): LOCATE P%+1%,1: PRINT TI$ 
NEXT 


TEMP2=1 / ( YMAX-YMTN ) : TEMP3=1/ ( XMAX-XMTN ) 
IF MRK$="" THEN MRK$="MARK  " 
REP2$=LEFTS (MRK$,1) 
T%=L%: REPS-LEFTS (CURVES(L$),1) 
IF T%>3 THEN T%=T%-3 : GOTO 1450 
IF REP$="C" THEN GOSUB 1870 
COUNT%=0: NDP%=0 
WHILE (PTR%<=NOD%) AND (X(PTR%) <= XMAX) 
IF X(PTR%) < XMIN THEN GOTO 1650 
OLDY%=Y%: OLDX%=x% 
NDP% = NDP%+1 
TEMP 1=ABS ( (X (PTR%) =XMIN) *TEMP3) *550 
TEMP=ABS ( (Y (PTR$)-YMIN) * TEMP2) *150 
DELTAY%=CINT (TEMP) : DELTAX%=CINT(TEMP1) 
Y%=TEMP4%-DELTAYS : X%=ORIGINX%+DELTAX3 
IF COUNT%=0 THEN OLDY%=Y%: OLDX%=X%: PSET(X%,Y%): 
FIRSTX%=X%: FIRSTY%=Y%: COUNT%= 
- IF REPS="P" THEN 1600 
LINE-(X%,Y%),1l 
IF REP2$ <> "M" AND REP <> "P" THEN GOTO 1630 
IF T$=1 THEN CIRCLE (X%,Y%),2,1,0,6.28,5/12 
IF T%=2 THEN LINE(X%,Y%)-(X%+6,Y%): LINE-(X%+3,V% 


-3): LINE-(X%,Y%) 

ELSE LINE(X%,Y%)-(X%+6,Y%-3),1,B : PSET(X%,Y%) 
GOTO 1640 

PSET(X%,Y%) 


IF REPS="C" THEN GOSUB 2140 
PTR3=PTR3+1 
WEND 

RETURN 


DIM H(200),A(200),L(200),U(200),Z(200) 

FOR 1% = 1 TO NOD%-1 

H(1%3)=X(13+1)-X(13) 

NEXT 

I%=2 TO NODS-1 

Tl=Y(I%+1)*H(I%-1) 

T2=Y (1%) *(X(I%+1)-X(I%-1)) 

T3=Y(I%-1)*H(I%) 

A(1%)=3*(T1-T2+T3)/(H(13-1)*H(13)) 

NEXT 13 

L(1)=1: U(l1)=0: 

FOR I%=2 TO NOD% 
L(13)=2*(X(13+1) -X(13-1)) -H(I3-1) *U(IS-1) 


FOR 


Z (1) =0 
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U(I%)=H(I% 
Z(I%)=(A(I 
NEXT I% 

L(NOD%)=1: Z(NOD%)=0: C(NOD%)=0 

FOR J%=NOD%-1 TO 1 STEP -1 

C(JT%)=Z (T%) -U(T%) *C(TS+1) 
T1=(Y(J$+1)-Y(J%))/H(J$) 
T2=H (3%) *(C(I%+1)+2*C(3%))/3 
B(J%)=T1-T2 
D(73)-(C(7X381)-C(73))/(3“H(72)) 
NEXT J% 

ERASE H,A,L,U,Z 


)/1(25) 
る) <H(I%-1) *Z(I%-1) ) /L(I%) 


IF PTR%+1 > NOD% OR X(PTR%+1) > XMAX THEN GOTO 2280 
T1=(XMAX-XMIN)/560 : T2=X(PTR%) 

TEMP1=ABS ( (X(PTR%+1) -XMIN) *TEMP3)*550 
DELTAX%=CINT ( TEMP1 ) 

XNEXT%=ORIGINX4+DELTAX% 

X%=X%+2 

IF X% >= XNEXTS 
T2=T2 + 2*T1 
T4=T2-X(PTR%) 
T3=Y(PTR%)+B(PTR%) *T4+0C (PTR%) *T4*T4+D(PTR%) *T4*T4*T4 
TEMP=ABS ( (T3-YMIN) *TEMP2) *150 

DELTAY%=CINT (TEMP): Y%=TEMP4%-DELTAY% 


THEN GOTO 2280 


IF Y%>172 THEN Y%=172 ELSE IF Y%<12 THEN Y%=12 
LINE-(X%,Y%) : GOTO 2190 

RETURN 

! ーーーーーー PLOT VITH GRIDS---------------------- 
FOR 1% = 1 TO 9 


T1%=ORIGINY$+I%*15 

LINE (73,T1%)-(623,T1%),,,6H4444. 
NEXT I% 

I$= 1 TO 13 

TIŻ-ORIGINXŻHIŻX40 
LINE(T1$,12)-(T13,162),,,SHAAAA 
NEXT I% 

RETURN 

END 


FOR 
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"NAME: Data Acquisition And Control (DAAC) 
| HEADER for BASICA 


"FILE NAME: DACHDR.BAS 
' 
.. DEVICE NAME: DAAC 


"RESERVED FUNCTION NAMES: 
AINM, AINS, AINSC, AOUM, AOUS, 
! BINM, BINS, BITINS, BITOUS, BOUM, BOUS, 
' CINM, CINS, CSET, DELAY 
"RESERVED DEF SEG VALUE NAME: DSEG 


(NAMES DEFINED AND USED BY HEADER: 
ADAPT%, AI, COUNT, FOUND% 
' HNAMES, sos, STATS 


'When using the BASICA Interpreter, this header 

'must be executed before any function calls are 

'made that access the DAAC adapter. It initializes 

'a number of variables for each function call. These 
'variables are reserved and should not be used except 
'to access the DAAC adapter. This routine also does a 
'DEF SEG to the segment where the DAAC Device Driver 
'(DAC.COM) is loaded. If you execute a DEF SEG to 
"access other hardware, you must DEF SEG to the segment 


'of the DAAC Device Driver before any subsequent 
‘calls to access the DAAC adapter. 

' 

' 

FOUND% = O 

SG% = &H2E 


'Start searching the interrupt vectors until you find 
'one that points to the DAAC device driver. 
"Do a DEF SEG to that segment. 
WHILE ((SG% <= &H3E) AND (FOUND% = 0)) 
DEF SEG = 0 
DSEG = PEEK(SG%) + PEEK(SG% + 1) * 256 
DEF SEG = DSEG 
HNAMES=" " 
FOR AI=10 TO 17 
HNAMES 2 HNAMES 4 CHRS(PEEK(AI)) 
NEXT AI 
IF HNAMEŞ = "DAAC " AND PEEK(18) + PEEK(19) <> 
O THEN FOUND% = 1 
SG% = SG% + 4 
WEND 
IF FOUND% = 0 THEN PRINT "ERROR: DEVICE DRIVER DAC.COM 
NOT FOUND" : END 
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"Nov initialize all function name variables for calls 
"to access the device driver. 

AINM = PEEK(&H13) * 256 + PEEK(8H12) 

AINS = PEEK(&H15) * 256 + PEEK(&H14) 

AINSC = PEEK(&H17) * 256 + PEEK(&H16) 

AOUM = PEEK(&H19) * 256 + PEEK(&H18) 

AOUS = PEEK(&H1B) * 256 + PEEK(&H1A) 

BINM = PEEK(SHID) * 256 + PEEK(SHLC) 

BINS = PEEK(&H1F) * 256 + PEEK(SHLE). 

BITINS = PEEK(8H21) * 256 + PEEK(8H20) 

BITOUS = PEEK(&H23) * 256 + PEEK(&H22) 

BOUM = PEEK(&H25) * 256 + PEEK(&H24) 

BOUS = PEEK(&H27) * 256 + PEEK(&H26) 

CINM = PEEK(&H29) * 256 + PEEK(&H28) 

CINS = PEEK(&H2B) * 256 + PEEK(&H2A) 

CSET = PEEK(&H2D) * 256 + PEEK(&H2C) 

DELAY = PEEK(&H2F) * 256 + PEEK(&H2E) 
'Finally, execute any call to re-initialize the 
'device driver from any former invocation of BASIC. 
ADAPT% = 0 


COUNT = 1 

STATS = 0 

CALL DELAY (ADAPT%, COUNT, STATS) 
' 


"End of DAAC BASICA Header 
t 


REM------ BALCAL.BAS (CALIBRATE BALANCE)---- 

REM 

DIM L(100) ,D(100) ,Y(100) ,PM(100) ,YM(100) ,RM(100) , LBS( 
100) 
DIM DAT(399) ,DAT%(399) ,DAT1(399) ,DAT1%(399) 


FOR TE=RISTO 12 
IF I = 1 THEN FILES = "LIFTP" 
IF I = 2 THEN FILES = "DRAGP" 
IF I = 3 THEN FILES = "YAWP" 
IF I = 4 THEN FILES = "PITCHMP" 
IF I = 5 THEN FILES = "YAWMP" 
IF I = 6 THEN FILES — "ROLLMP" 
IF I = 7 THEN FILES = "LIFTN" 
IF I = 8 THEN FILES = "DRAGN" 
IF I = 9 THEN FILES = "YAWN" 
IF I = 10 THEN FILES = "PITCHMN" 
IF I = 11 THEN FILES = "YAWMN" 
IF I = 12 THEN FILES = "ROLLMN" 
REM 
REM 
GOSUB 1380 ‘RECORD CALIBRATION DATA 
NEXT I 
CES 


LOCATE 10,5: PRINT "CALCULATING CALIBRATION CONSTANTS, 
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PLEASE VAIT" 


960 GOSUB 2100 "CALCULATE KISK2 FOR ALL PRIME GAGES 

970 GOSUB 2640 'CONVERT COUNTS TO FORCES 

980 GOSUB 3180 'CALCULATE INTERACTION CONSTANTS 

990 REM 

1000 REM STORE CALIBRATION CONSTANTS 

1010 OPEN "C:CONST" FOR OUTPUT AS #1 

1020 WRITE 41, INCALL, INCALD, INCALY , INCALPM, INCALYM, INCALRM 

1040 WRITE fl,KILPOS,K2LPOS,KIDPOS,K2DPOS,KIVPOS,K2VPOS 

1050 WRITE #1,K1PMPOS,K2PMPOS ,K1YMPOS,K2YMPOS,K1RMPOS, 
K2RMPOS | 

1060 WRITE #1,K1LNEG, K2LNEG,K1DNEG,K2DNEG, K1YNEG, K2YNEG 

1070 WRITE f1,KIPMNEG,K2PMNEG,KIYMNEG,K2YMNEG,KIRMNEG, 
K2RMNEG | 

1080 WRITE +1,DDDL1P,DDDL2P,DYDL1P,DYDL2P,DPMDL1P,DPMDL2P 
,DYMDL1P, DYMDL2P 

1090 WRITE #1,DRMDL1P,DRMDL2P 

1100 WRITE #1,DLDD1P,DLDD2P,DYDD1P,DYDD2P,DPMDD1P,DPMDD2P 
,DYMDD1P, DYMDD2P 

1110 WRITE #1,DRMDD1P,DRMDD2P 

1120 WRITE $Ż1,DLDY1P,DLDY2P,DDDY1P,DDDY2P,DPMDY1P,DPMDY2P 
, DYMDY1P,DYMDY2P 

1130 WRITE #1,DRMDY1P,DRMDY2P 

1140 WRITE #1,DLDPM1P,DLDPM2P,DDDPM1P,DDDPM2P,DYDPM1P,DYD 
PM2P,DYMDPMIP,DYMDPM2P 

1150 WRITE f1,DRMDPMIP,DRMDPM2P 

1160 WRITE #1,DLDYM1P, DLDYM2P,DDDYM1P, DDDYM2P, DYDYM1P, DYD 
YM2P, DPMDYM1P, DPMDYM2P 

1170 WRITE #1,DRMDYM1P,DRMDYM2P 

1180 WRITE #1,DLDRM1P,DLDRM2P,DDDRM1P,DDDRM2P,DYDRM1P,DYD 
RM2P,DPMDRMIP,DPMDRM2P 

1190 WRITE #1,DYMDRM1P, DYMDRM2P 

1200 WRITE $41,DDDL1N,DDDL2N,DYDL1N,DYDL2N,DPMDL1N,DPMDL2N 
, DYMDL1N , DYMDL2N 

1210 WRITE $41,DRMDL1N, DRMDL2N 

1220 WRITE #1,DLDD1N, DLDD2N,DYDD1N,DYDD2N,DPMDD1N, DPMDD2N 
, DYMDD1N , DYMDD2N 

1230 WRITE #1,DRMDD1N, DRMDD2N | 

1240 WRITE #1,DLDY1N,DLDY2N,DDDY1N,DDDY2N,DPMDY1N,DPMDY2N 
, DYMDY1N , DYMDY2N 

1250 WRITE $41,DRMDY1N,DRMDY2N 

1260 WRITE f1,DLDPMIN,DLDPM2N,DDDPMIN,DDDPM2N,DYDPMIN,DYD 
PM2N,DYMDPMIN 

1270 WRITE #1,DYMDPM2N,DRMDPM1N, DRMDPM2N 

1280 WRITE #1,DLDYM1N,DLDYM2N,DDDYM1N,DDDYM2N,DYDYM1N,DYD 
YM2N , DPMDYM1N 

1290 WRITE #1,DPMDYM2N,DRMDYM1N, DRMDYM2N 

1300 WRITE +41, DLDRM1N, DLDRM2N, DDDRM1N, DDDRM2N, DYDRM1N, DYD 
RM2N , DPMDRM1N 

1310 WRITE #1,DPMDRM2N,DYMDRM1N,DYMDRM2N 
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CLOSE #1 


CLS: LOCATE 10,5 

PRINT " CALIBRATION COMPLETE" 
END 

REM 

REM 


REM---RECORD CALIBRATION DATA----== 

COLOR 15,1: tr CTS 

PRINT "CALIBRATION FOR THE LOADING OF THE ";FILES;" 
COMPONENT" 


PRINT: INPUT "AFTER ALL AMPLIFIERS ARE ZEROED PRESS 
RETURN" ;X 

CLS 

PRN 21 ZL ZY ZPM ZYM 

ZRM " 

STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 


CALL AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 
, STOR% , COUNT, RATE, DAT% (0) , STATS) | 

ZD = 0:Z2PM=0: ZL=0: ZYM=0 

IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

FOR J = 0 TO 396 STEP 4 

DAT (J) =(DAT% (J) /204.8) -10 

ZD =ZD + DAT(J) 


NEXT J 
ZD = ZD/100 
FOR J = 1 TO 397 STEP 4 


DAT (J)=(DAT% (J)/204.8)-10 


ZES, ZE Ји 
NEXT J 
FOR J = 2 TO 398 STEP 4 


DAT (J) =(DAT% (J) /204.8)-10 


ZPM = ZPM + DAT(J) 
NEXT J 
FOR J = 3 TO 399 STEP 4 


DAT (J)=(DAT% (J)/204.8)-10 


ZYM = ZYM + DAT(J) 

NEXT J 

STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT$= 1: DEVICES- CHANLO%= 0: CTRL%= 0: CHANHI%=1 


CALL AINSC(ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 
, STOR る , COUNT, RATE, DAT1% (0) , STATS) 

ZY =0:ZRM=0 

IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STAT%:END 

FOR J = 0 TO 198 STEP 2 

DAT1 (J) =(DAT1% (J) /204.8)-10 

ZRM = ZRM + РАТ1(Ј) 

NEXT J 
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ZRM = ZRM/100 

FOR J = 1 TO 199 STEP 2 
DAT1(J)=(DAT1%(J)/204.8)-10 
ZY = ZY + DAT1(J) 

NEXT J 


ZL=ZL/100:ZPM=ZPM/100:ZYM=ZYM/100:ZY=ZY/100 

LOCATE 3,1: PRINT USING "+#.###":Z2D: LOCATE 3,10: PRI 
NT USING fl-4, £4HĦ ZI, 

LOCATE 3,19: PRINT USING "+#.###";ZY 

LOCATE 3,28: PRINT USING "+#.###";ZPM 

LOCATE 3,37: PRINT USING "+#.###";ZYM 

LOCATE 3,46: PRINT USING "+#.###";ZRM 

REM 

PRINT: PRINT 

INPUT "AFTER PLACING ALL CAL SWITCHES TO + SETTING HIT 
RETURN" ; X 

PRINT 

PRINT " CALD CEL CALY CALPM CAŁYM 
CALRM" 

Su = MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 0: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 
CALL AINSC(ADAPT%, DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 


, STOR% , COUNT , RATE, DAT% (0) , STATS) 
CALD=0 : CLL=0 : CALYM=0 : CALPM=0 

IF STAT%<> 0,THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 

FOR J = 0 TO 396 STEP 4 

DAT (J) =(DAT% (J) /204.8)-10 

CALD =CALD + DAT(J) 

NEXT J 

CALD = CALD/100 

FOR J = 1 TO 397 STEP 4 

DAT (J) =(DAT% (J) /204.8)-10 

CLL = CLL + DAT(J) 

NEXT J 

FOR J = 2 TO 398 STEP 4 

DAT (J) =(DAT% (J) /204.8)-10 


CALPM = CALPM + DAT(J) 

NEXT J 

FOR J = 3 TO 399 STEP 4 

DAT (J) =(DAT% (J) /204.8)=-10 

CALYM = CALYM + DAT(J) 

NEXT J 

STATŻ= MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 1: DEVICE%= 9: CHANLOŻ= 0: CTRL$= 0: CHANHI%= 
CALL AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL$%,MODE% 


, STOR%, COUNT, RATE, DAT1% (0) ,STAT$) 

CALY =0:CALRM=0 

IF STAT%<> 0 THEN PRINT USING "EXECUTION ERROR ###"; 
STATS: END 
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FOR J = 0 TO 198 STEP 2 
DAT1 (J) =(DAT1% (J) /204.8)-10 
CALRM = CALRM + DAT1(J) 
NEXT J 

CALRM = CALRM/100 

FOR J = 1 TO 199 STEP 2 
DAT1(J)=(DAT1%(J)/204.8)-10 
CALY = CALY + DAT1(J) 


NEXT J 
CLL=CLL/100:CALPM=CALPM/100:CALYM=CALYM/100:CALY=CALY 
/100 
LOCATE 10,1: PRINT USING "+#.###";CALD 
LOCATE 10,10: PRINT USING "+#.###";CLL 
LOCATE 10,19: PRINT USING "+#.###">;CALY 
LOCATE 10,28:. PRINT USING "+#.###";CALPM 
LOCATE 10,37: PRINT USING "+#.###";CALYM 
LOCATE 10,46: PRINT USING "+#.###";CALRM 
REM 
REM 

IF I = 1 THEN INCALL = CLL- Z 

IF I = 2 THEN INCALD = CALD 7728 

IF I = 3 THEN INCALY = CAIT 2777 

IF I = 4 THEN INCALPM = CALPM - ZPM 

IF I = 5 THEN INCALYM = CALYM - ZYM 

IF I = 6 THEN INCALRM = CALRM - ZRM 


PRINT: INPUT "PRESS ENTER (RETURN) TO CONTINUE";X 
CLS:PRINT "RETURN ALL CAL SWITCHES TO CENTER POSITION" 
PRINT: | 

PRINT "LOAD THE ";FILES;" GAGE FROM O TO 20 POUNDS AND 
THEN BACK TO 0" 

PRINT "POUNDS IN 1 POUND INCREMENTS. 
TWO READINGS FOR" | 
PRINT "EACH POUND WEIGHT EXCEPT FOR THE LOAD AT 20 
POUNDS WHICH WILL HAVE" 


THERE SHOULD BE 


PRINT "ONLY ONE READING" 

PRINT 

PRINT "PRESS THE Fl KEY WHEN THE LOADING IS FINISHED" 
PRINT "PRESS THE F2 KEY WHEN READY TO RECORD THE DATA 


FOR THAT LOAD" 


PRINT: INPUT "PRESS ENTER (RETURN) TO CONTINUEN, X 


CLS 

PRINT " DRAG LIFT YAW PITCH YAW 
ROLL " 

PRINT " MOM. MOM. 
MOM. " 

SOAP = 0: N=4 

FOR K = 1 TO 100 

ON KEY(1) GOSUB 1830 (SET STOP FLAG 

ON KEY(2) GOSUB 1841 (RECORD DATA 

KEY(1) ON: KEY(2) ON 
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IF SOAP = 2 THEN GOTO 1780 

IF SOAP = 1 THEN GOTO 1920 

GOTO 1730 

SOAP = 0 

NOD% = K 

NEXT K 

GOTO 1920 

REM SET STOP FLAG 

SOAP = 1 

RETURN 

REM STEPS TO RECORD DATA 

STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 0: DEVICES= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=3 
CALL AINSC (ADAPT%, DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 


, STOR% , COUNT , RATE , DAT% (0) , STATS) 
D(K)=0:L(K)=0:YM(K)=0: PM(K) =0 

IF STAT%<> O,THEN PRINT USING "EXECUTION ERROR ###"; 
STAT% : END 


FOR J = 0 TO 396 STEP 4 
DAT (J) =(DAT% (J) /204.8)=-10 
D(K) =D(K) + DAT(J) 

NEXT J 

D(K) = D(K)/100 

FOR J = 1 TO 397 STEP 4 
DAT(J)=(DAT%(J)/204.8)-10 


L(K) = L(K) + DAT(J) 
NEXT J 

FOR J = 2 TO 398 STEP 4 
DAT (J)=(DAT$ (J)/204.8)-10 


PM(K) = PM(K) + DAT(J) 

NEXT J 

FOR J = 3 TO 399 STEP 4 

DAT (J) =(DAT% (J) /204.8)-10 

YM(K) = YM(K) + DAT(J) 

NEXT J 

STAT%=0: MODE%=0: STOR%=0: COUNT=100: RATE=500 
ADAPT%= 1: DEVICE%= 9: CHANLO%= 0: CTRL%= 0: CHANHI%=1 
CALL AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL%,MODES 


, STORS , COUNT, RATE, DAT1% (0) , STATS) 

(K) =0:RM(K)=0 

IF STAT%<> O THEN PRINT USING "EXECUTION ERROR ###"; 
STATS: END 


FOR J = 0 TO 198 STEP 2 
DAT1(J)=(DAT1%(J)/204.8)-10 


RM(K) = RM(K) + DAT1(J) 
NEXT J 

RM(K) = RM(K)/100 

FOR J = 1 TO 199 STEP 2 
DAT1(J)=(DAT1%(J)/204.8)-10 


Y(K) = Y(K) + DATL(J) 
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NEXT J 
L(K)=L(K)/100:PM(K)=PM (K) /100:YM (K) =YM (K) /100:Y (K)=Y ( 
K)/100 
IF N> 23 THEN N=4:CLS:PRINT " DRAG LIFT YAW 
PITCH YAW ROLL ":PRINT " 
MOM. MOM. MoM. " 
LOCATE N,1: PRINT USING "+#.###";D(K) 
LOCATE N,10: PRINT USING "+#.###";L(K) 
LOCATE N,19: PRINT USING "+#.###";Y(K) 
LOCATE N,28: PRINT USING "+#.###";PM(K) 
LOCATE N,37: PRINT USING "+#.###";YM(K) 
LOCATE N,46: PRINT USING "+#.###";RM(K) 
N=N+1 
SOAP = 2 
RETURN 'GO BACK FOR NEXT DATA POINT 


"---CORRECT FOR DRIFT------------- 
DIFFL =(L(NOD%)-L(1)): DIFFD=(D(NOD%)-D(1)):DIFFY=(Y( 
NOD%) -Y(1)) 
DIFFPM=(PM(NOD%) -PM(1)) : DIFFYM=(YM(NOD%) -YM(1)) :DIFFR 
M= (RM (NOD%) -RM(1) ) 
D=NOD%-1 
CORL=DIFFL/D:CORD=DIFFD/D:CORY=DIFFY/D:CORPM=DIFFPM/D 
: CORYM=DIFFYM/D:CORRM=DIFFRM/D 
Azı 
FOR K= 2 TO NOD% 

L(K)=L(K) - (A*CORL) 

D(K) =D(K) - (A*CORD) 

Y (K)=Y(K)-(A*CORY) 

PM (K) =PM(K) - (A*CORPM) 

YM (K) =YM (K) - (A*CORYM) 

RM (K) =RM (K) - (A*CORRM) 

A = A+1 
NEXT K | 
"---VRITE DATA TO FILE----- 
FILES = "C:"+FILES 
OPEN FILES FOR OUTPUT AS #1 
VRITE 41,ZL,ZD,ZY,ZPM,ZYM,ZRM,NOD$ 
VRITE 41,CLL,CALD,CALY,CALPM,CALYM,CALRM 


IF I=1 OR I=2 OR I=3 THEN GOTO 2079 
IF I=7 OR I=8 OR I=9 THEN GOTO 2067 
LBS = 
ORM TE AI 

IF I=4 OR I=10 THEN D1=20 

IF I=5 OR I=ll THEN D1=4.5 

IF I=6 OR I=12 THEN D1=11.5 

MOM = LBS*D1 'CONVERT TO MOMENTS 


IF I=10 OR I=11 OR I=12 THEN MOM=MOM* (-1) 
WRITE #1,L(J),D(J),Y(J),PM(J) , YM(J) ,RM(J) , MOM 
LBS = LBS+1 

NEXT J 
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I 


LBS 
NEXT J 
CLOSE #1 


J 
PE 
E 


RETURN 


LBS = 


FOR J 


0 


— 
— 


2 TO NOD% 


OR I=10 THEN D1=20 
OR I=11 THEN D1=4.5 
IF I=6 OR I=12 THEN D1=11.5 
MOM 
IF I=10 OR I=11 OR I=12 THEN MOM=MOM* (-1) 
WRITE $1,L(J),D(J),Y(J),PM(J),YM(J) , RM(J) , MOM 


LBS*D1 


LBS=1 


! CONVERT TO MOMENTS 


!GO BACK TO SET UP NEXT CALIBRATION 


1 TO 21 
WRITE $1,L(J),D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS 


JBS = DBS = 1 
NEAT J 

LBS = =19 

FOR J = 22 TO NOD% 


WRITE #1,L(J),D(J),Y(J),PM(J),YM(J),RM(J),LBS 


LBS 


NEXT J 


CLOSE #1 
RETURN 


LBS = 
FOR J 


0 


LBS 1 


"GO BACK TO SET UP NEXT CALIBRATION 


1 TO 21 
WRITE #1,L(J),D(J),¥(J),PM(J) , YM(J) ,RM(J), LBS 


LBS EBS + h 
NEXT J 

LBS = 19 

FOR J = 22 TO NOD% 


WRITE #1,L(J),D(J),Y(J),PM(J) ,YM(J) ,RM(J) , LBS 


LBS LBS - 1 

NEXT J 
CLOSE +1 
RETURN 'GO BACK TO SET UP NEXT CALIBRATION 
REM--------------------------------------------- 
REM CALCULATE K1&K2 FOR PRIME GAGES 
FOR I = 1 TO 12 

DB — 1 THEN FLILEŞ = "LIFTP" 

IF I = 2 THEN FILE$ = "DRAGP" 

IF I — 3 THEN FILEŞ — "YAVP" 

IF I = 4 THEN FILES = "PITCHMP" 

IF I = 5 THEN FILES = "YAWMP" 

IF I = 6 THEN FILES = "ROLLMP" 

IF I = 7 THEN FILE$ = "LIFTN" 

IF I = 8 THEN FILE$ = "DRAGN" 

IF I = 9 THEN FILES = "YAWN" 

IF I = 10 THEN FILES = "PITCHMN" 

IF I = 11 THEN FILES = "YAWMN" 


Figure A.16 
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IF I = 12 THEN FILE$ = "ROLLMN" 
GOSUB 2295 "READ FILES AND PERFORM CALCULATIONS 
NEXT L 
RETURN IGO BACK TO CONVERT COUNTS 
REM 
REM 


A=0 : B=0: C=0: D=0: E=0: 

OPEN "C:"+FILES FOR INPUT AS #2 

INPUT %2,ZL,ZD,ZY,ZPM,ZYM, ZRM,NOD% 

INPUT #2,CLL,CALD, CALY, CALPM, CALYM, CALRM 


FOR J = 1 TO NOD% 
INPUT 42,L(J),D(J),Y(J) ,PM(J),YM(J) ,RM(J) ,LBS (J) 
IF I = 1 OR I = 7 THEN X = L(J) 
IF I = 2 OR I = 8 THEN X = D(J) 
IF I = 3 OR I = 9 THEN X = Y(J) 
IF I = 4 OR I = 10 THEN X = PM(J) 
IF I = 5 OR I = 11 THEN X = YM(J) 
IF I = 6 ORI = 12 THEN X = RM(J) 
А = A + (X42) 
B = B + (XA3) 
C =C + (X^ 
D = D + (X*LBS(J)) 
E = E + (X*X*LBS(J)) 
NEXT J 
CLOSE #2 
K2 = ((DZ/A)-(E/B)) / ( (B/A)—-(C/B) ) 
Kl = (D/A) - (K2*(B/A)) 
IF I = 1 THEN K1LPOS = Kl: K2LPOS = K2 
IF I = 2 THEN K1DPOS = Kl: K2DPOS = K2 
IF I = 3 THEN K1YPOS = Kl: K2YPOS = K2 
IF I = 4 THEN K1PMPOS = Kl: K2PMPOS = K2 
IF I = 5 THEN K1YMPOS = Kl: K2YMPOS = K2 
IF I = 6 THEN K1RMPOS = Kl: K2RMPOS = K2 
IF I = 7 THEN K1LNEG = Kl: K2LNEG = K2 
IF I = 8 THEN KIDNEG = Kl: K2DNEG = K2 
IF I = 9 THEN KIYNEG = Kl: K2YNEG = K2 
IF I = 10 THEN K1PMNEG = Kl: K2PMNEG = K2 
IF I = 11 THEN K1YMNEG = Kl: K2YMNEG = K2 
IF I = 12 THEN K1RMNEG = Kl: K2RMNEG = K2 
RETURN  'GO BACK TO CALCULATE K1&K2 FOR NEXT FILE 
REM----------------------------------------------- 
REM CONVERT COUNTS TO FORCES 
FOR I = 1 TO 12 
IF I = 1 THEN FILES = "LIFTP" 
IF I = 2 THEN FILES = "DRAGP" 
IF I = 3 THEN FILES = "YAWP" 
IF I = 4 THEN FILES = "PITCHMP" 
IF I = 5 THEN FILES = "YAWMP" 
IF I = 6 THEN FILES = "ROLLMP" 
IF I = 7 THEN FILES = "LIFTN" 
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IF I — 8 THEN FILEŞ — "DRAGN" 
IF I = 9 TREN FILES = "YAWN" 
IF I = 10 THEN FILES = "PITCHMN" 
IF I = 11 THEN FILES = "YAWMN" 
Tr = 12 TREN FILES — "ROLLU" 
GOSUB 2835 'READ IN FILE AND CONVERT 
NEXT I 
RETURN 'GO BACK TO CALCULATE INTERACTON CONSTANTS 
REM 
REM 


A=0 : B=0: C=0: D=0: E=0: F=0 

OPEN "C:"+FILES FOR INPUT AS #2 

INPUT +2,ZL,ZD,ZY,ZPM,ZYM, ZRM, NOD% 

INPUT 42,CLL,CALD,CALY,CALPM, CALYM, CALRM 

FOR J = 1 TO NOD% 
INPUT 42,L(J) ,D(J) ,Y(J),PM(J) ,YM(J) ,RM(J) ,LBS(J) 
A=((INCALL/CLL-ZL) *(L(J)-ZL)) 
B= ( (TNCALD/CALD-ZD) *(D( け ) -ZD) ) 
C-( (INCALY/CALY-ZY) “ (Y(7)-ZY)) 
D=((INCALPM/CALPM-ZPM) *(PM(J)-ZPM)) 
E=((INCALYM/CALYM-ZYM) *(YM(J)-ZYM)) 
F=((INCALRM/CALRM-ZRM) * (RM(J)-ZRM) ) 
IF L(J) < 0 THEN Kl = K1LNEG: K2 = K2LNEG ELSE Kl 
= KILPOS:K2 = K2LPOS 
L(J)=(K1*A) + (K2*(A72)) 


IF D(J) < 0 THEN Kl = K1DNEG: K2 = K2DNEG ELSE Kl 
= KIDPOS:K2 = K2DPOS 

D(J)=(K1*B) + (K2*(B"2)) 

IF V(J) < 0 THEN Kl = K1YNEG: K2 = K2YNEG ELSE Kl 


= K1YPOS:K2 = K2YPOS 
Y(J)=(K1*C) + (K2*(C0%2)) 
IF PM(J) < O THEN Kl=K1PMNEG:K2=K2PMNEG ELSE K1=K1 
PMPOS : K2=K2PMPOS 
PM(J)=(K1*D) + (K2*(D"2)) 
IF VM(J) < O THEN Kl=K1YMNEG:K2=K2YMNEG ELSE Kl=K1 
YMPOS : K2=K2YMPOS 
YM(J)=(K1*E) + (K2*(E"2)) 
IF RM(J) < 0 THEN Kl=K1RMNEG:K2=K2RMNEG ELSE K1=K1 
RMPOS : K2=K2RMPOS 
RM(J)=(K1*F) + (K2*(F*2)) 
NEXT J 
CLOSE #2 
REM SAVE CONVERTED COUNTS 
OPEN "B:"+FILES FOR OUTPUT AS #1 
WRITE 41,ZL,ZD,ZV,ZPM,ZVM,ZRM,NODŻ 
WRITE #1,CLL,CALD, CALY, CALPM, CALYM, CALRM 


FOR J = 1 TO NOD% 
WRITE 41,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS(J) 
NEXT J 

CLOSE #1 
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RETURN "GO BACK TO CALCULATE INTERACTION CONSTANTS 


FOR I = 1 16” 

IF I = 1 THEN FILES = "LIFTE" 

IF I = 2 THEN FILES = "DRAGP" 

IF I = 3 THEN FILES = "YAWP" 

IF I = 4 THEN FILES — "PITCHMP" 

IF I = 5 THEN FILES = "YAWMP" 

IF I = 6 THEN FILES = "ROLLMP" 

IF I = 7 THEN FILES = "LIFTN" 

IF I = 8 THEN FILES = "DRAGN" 

IF I = 9 THEN FILES = "YAWN" 

IF I = 10 THEN FILES = "PITCHMN" 

IF I = 11 THEN FILES = "YAWMN" 

IF I = 12 THEN FILES = "ROLLMN" 
GOSUB 3365 ‘READ IN FILES AND CALCULATE CONSTANTS 
NEXT 1 
RETURN "GO BACK TO SAVE CALIBRATION DATA 
REM 
REM 


A=0: B=0: C=0: D1=0: D2=0: D3=0: D4=0: D5=0: E1=0: E2=0: E3=0:E 
4=0:E5=0 
OPEN "B:"+FILES FOR INPUT AS #2 
INPUT ;Ż2,ZL,ZD,ZY,ZPM,ZYM,ZRM,NOD% 
INPUT 42,CLL,CALD,CALV,CALPM,CALVM,CALRM 
FOR J = 1 TO NOD% 

INPUT 2,L(J),D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS(J) 

IF I=1 OR I=7 THEN X=L(J):Y1=D(J):Y2=Y(J):Y3=PM(J) 
:Y4=YM(J) :Y5=RM(J) 

IF I=2 OR I=8 THEN X=D(J):Y1=L(J):Y2=Y(J):Y3=PM(J) 
: Y4=YM(J) : Y5=RM(J) 

IF I=3 OR I=9 THEN X=Y(J):Y1=L(J):Y2=D(J):Y3=PM(J) 
:Y4=YM(J) :Y5=RM(J) | 

IF I=4 OR I=10 THEN X=PM(J):Y1=L(J):Y2=D(J) :Y3=Y(J) 
: Y4=YM(J) : Y5=RM(J) 

IF I=5 OR I=11 THEN X=YM(J):Y1=L(J):Y2=D(J):Y3=Y(J) 
:Y4=PM(J) :Y5=RM(J) . 

IF I=6 OR I=12 THEN X=RM(J):Y1=L(J):Y2=D(J):Y3=Y(J) 
:Y4=PM(J) :YS=YM(J) 

A € A 4 (X”2) 


B = B  (X”3) 

C = C + (X^4) 

О1 = О1 + (X*Y1) 

El = El + ((X"2)*Y1) 
D2 = D2 + (X*Y2) 

Е2 — E2 1 ((X”2)Y2) 
D3 = D3 + (X*Y3) 

E3 = E3 + ((X"2)*Y3) 
D4 = D4 + (XxY4) 
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E4 = E4 + ((X^2)*Y4) 

D5 = D5 + (X*Y5) 

E5 = E5 + ((X^2)*Y5) 

NEXT J 
CLOSE #2 
KE DE/A)-(EL/B))/ ((8. R) 
К11 = (О1/А)-(К12*(В/А)) 
177” ((D2/Z/A) S (E2ZB))/ ( (BZAÀ) PB 
K21 = (D2/A)-(K22*(B/A)) 
5775 HU) ” (E375) )/ ( (B87 A) =( C78) ) 
K31 = (D3/A)-(K32*(B/A)) 
K42 = ((D4/A)-(E4/B))/((B/A)-(C/B)) 
K41 = (D4/A)-(K42*(B/A)) 
Eu (D5/A)=(ES/B))/ ((B/A)—-(C/B) ) 
K51 = (D5/A)-(K52*(B/A)) 


IF I=1 THEN DDDL1P=K11:DDDL2P=K12:DYDL1P=K21:DYDL2P=K 
22 :DPMDL1P=K31 

IF I=1 THEN DPMDL2P=K32 : DYMDL1P=K41 : DYMDL2P=K4 2 : DRMDL 
1P=K51:DRMDL2P=K52 

IF I=2 THEN DLDD1P=K11: DLDD2P=K12 : DYDD1P=K21: DYDD2P=K 
22:DPMDD1P=K31 

IF I=2 THEN DPMDD2P=K32:DYMDD1P=K41:DYMDD2P=K42: DRMDD 
1P=K51: DRMDD2P=K52 

IF I=3 THEN DLDY1P=K11:DLDY2P=K12:DDDY1P=K21:DDDY2P=K 
22:DPMDY1P=K31 

IF 1=3 THEN DPMDY2P=K32: DYMDY 1P=K41:DYMDY2P=K4 2 : DRMDY 
1P=K51:DRMDY2P=K52 

IF I=4 THEN DLDPM1P=K11: DLDPM2P=K12 : DDDPM1P=K21: DDDPM 
2P=K22 : DYDPM1P=K3 1 

IF I=4 THEN DYDPM2P=K32 : DYMDPM1P=K41: DYMDPM2 P=K42 : DRM 
DPM1P=K51:DRMDPM2P=K52 

IF I=5 THEN DLDYM1P=K11: DLDYM2 P=K12 : DDDYM1P=K21:DDDYM 
2P=K22:DYDYMIP=K31 

IF 1-5 THEN DYDYM2PFK32:DPMDYMIP-K41:DPMDYM2P-K42:DRM 
DYM1P=K51: DRMDYM2 P=K52 

IF I=6 THEN DLDRM1P=K11: DLDRM2P=K12 : DDDRM1P=K21: DDDRM 
2P=K22:DYDRM1P=K3 1 

IF 1=6 THEN DYDRM2P=K3 2 : DPMDRM1P=K4 1: DPMDRM2P=K4 2 : DYM 
DRM1P=K51 : DYMDRM2 P=K52 

IF I=7 THEN DDDL1IN=K11: DDDL2N=K12 : DYDLIN=K21: DYDL2N=K 
22: DPMDLIN=K31 

IF I=7 THEN DPMDL2N=K32:DYMDL1N=K41:DYMDL2N=K42:DRMDL 
1N=K51:DRMDL2N=K52 

IF I=8 THEN DLDDIN=K11: DLDD2N=K12: DYDDIN=K21: DYDD2N=K 
22:DPMDDIN=K31 

IF I=8 THEN DPMDD2N=K3 2: DYMDDIN=K41: DYMDD2N=K42 : DRMDD 
1N=K51: DRMDD2N=K52 

IF 1-9 THEN DLDY1N=K11:DLDY2N=K12:DDDY1N=K2 1:DDDY2N=K 
22:DPMDY1N=K31 

IF I=9 THEN DPMDY2N=K32:DYMDY1N=K41:DYMDY2N=K42: DRMDY 
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1N=K51:DRMDY2N=K52 

IF I=10 THEN DLDPM1N=K11: DLDPM2N=K12 : DDDPM1N=K21:DDDP 
M2N=K22:DYDPM1N=K31 

IF I=510 THEN DYDPM2N=K32:DYMDPM1N=K41:DYMDPM2N=K42 :DR 
MDPM1N=K51 

IF I=10 THEN DRMDPM2N=K52 

IF I=11 THEN DLDYM1N=K11: DLDYM2N=K12 : DDDYM1IN=K21: DDDY 
M2N-K22:DYDYMIN-K31 

IF I=11 THEN DYDYM2N=K32 : DPMDYM1IN=K41: DPMDYM2N=K42:DR 
MDYM1N=K51 

IF I=11 THEN DRMDYM2N=K52 

IF I=12 THEN DLDRM1N=K11: DLDRM2N=K12 : DDDRM1IN=K21:DDDR 
M2N-K22:DYDRMIN-K31 

IF I=12 THEN DYDRM2N=K32 : DPMDRM1IN=K41: DPMDRM2N=K42: DY 
MDRM1N=K51 

IF I-12 THEN DYMDRM2N=K52 

RETURN "GOBACK TO SAVE CALIBRATION DATA 
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APPENDIX B 


TABLES 


SEES ET Ure OR TEE TQDEE 62559 TRANSDUCER AMPLIFIER 


Amp. # Component Gain Set. Filter Set. Exe Volt: 
l DRAG MAX VAR. ıl 55.0 
2 EE MAX VAR. b "6.5 
3 PTICH M. EŃ ii 0 
4 YAW M. MAX VAR. al 250 
5 ROLL M. MAX VAR. aL 55.0 
6 YAW MAX VAR. Jl +5 10 
Y AOA ır JU “50 


Table B.1 Amplifier Set-Up 
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